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1 Introduction 

QCD as a theory of extended, strongly interacting matter is familiar from big 
bang evolution which, within the time interval from electro-weak decoupling 
(10~^^ s) to hadron formation (5 • 10^® s), is dominated by the expansion 
of quark-gluon matter, a color conducting plasma that is deconfined. In the 
1970's deconfinement was predicted [1, 2, 3] to arise from the newly discovered 
"asymptotic freedom" property of QCD; i.e. the plasma was expected to be a 
solution of perturbative QCD at asymptotically high square-momentum trans- 
fer Q^, or temperature T. Thus the Quark-Gluon-Plasma (QGP) was seen as 
a dilute gas of weakly coupled partons. This picture may well hold true at 
temperatures in the GeV to TeV range. However it was also known since R. 
Hagedorns work [4] that hadronic matter features a phase boundary at a very 
much lower temperature, T{H) = 170 McV. As it was tempting to identify 
this temperature with that of the cosmological hadronization transition, thus 
suggesting T{H) = T(crit), the QGP state must extend downward to such a 
low temperature, with « IGeV^, and far into the non-perturbative sector 
of QCD, and very far from asymptotic freedom. The fact that, therefore, the 
confinement-deconfinement physics of QCD, occuring at the parton-hadron 
phase boundary, had to be explained in terms other than a dilute perturba- 
tive parton gas, was largely ignored until rather recently, when laboratory 
experiments concerning of the QGP had reached maturity. 

In order to recreate matter at the corresponding high energy density in 
the terrestial laboratory one collides heavy nuclei (also called "heavy ions" ) at 
ultrarelativistic energies. Quantum Chromo dynamics predicts [2, 5, 6] a phase 
transformation to occur between deconfined quarks and confined hadrons. At 
near-zero net baryon density (corresponding to big bang conditions) non- 
perturbative Lattice-QCD places this transition at an energy density of about 
1 GeV/fm^, and at a critical temperature, Tcrti » 170 MeV [6, 7, 8, 9, 10]. 
The ultimate goal of the physics with ultrarelativistic heavy ions is to locate 
this transition, elaborate its properties, and gain insight into the detailed 
nature of the deconfined QGP phase that should exist above. What is meant 
by the term " ultrarelativistic" is defined by the requirement that the reaction 
dynamics reaches or exceeds the critical density e « 1 GeV/ fw?. Required 
beam energies turn out [10] to be ^/s > 10 GeV , and various experimental 
programs have been carried out or are being prepared at the CERN SPS (up 
to about 20 GeV), at the BNL RHIC colUder (up to 200 GeV) and finally 
reaching up to 5.5 TeV at the LHC of CERN. 

QCD confinement-deconfinement is of course not limited to the domain 
that is relevant to cosmological expansion dynamics, at very small excess of 
baryon over anti-baryon number density and, thus, near zero baryo-chemical 
potential /j-b- In fact, modern QCD suggests [11, 12, 13] a detailed phase di- 
agram of QCD matter and its states, in the plane of T and baryo-chemical 
potential /is. For a map of the QCD matter phase diagram we are thus em- 
ploying the terminology of the grand canonical Gibbs ensemble that describes 
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an extended volume V of partonic or hadronic matter at temperature T. In 
it, total particle number is not conserved at relativistic energy, due to parti- 
cle production-annihilation processes occurring at the microscopic level. How- 
ever, the probability distributions (partition functions) describing the particle 
species abundances have to respect the presence of certain, to be conserved 
net quantum numbers (i), notably non-zero net baryon number and zero net 
strangeness and charm. Their global conservation is achieved by a thermody- 
namic trick, adding to the system Lagrangian a so-called Lagrange multiplier 
term, for each of such quantum number conservation tasks. This procedure 
enters a "chemical potential" /Xj that modifies the partition function via an 
extra term exp (—fii/T) occuring in the phase space integral (see section 4 for 
detail). It modifies the canonical "punishment factor" exp(— i?/T), where E 
is the total particle energy in vacuum, to arrive at an analogous grand canon- 
ical factor for the extended medium, of cxp (—E/T — fii/T). This concept is 
of prime importance for a description of the state of matter created in heavy 
ion collisions, where net-baryon number (valence quarks) carrying objects are 
considered extended "fireballs" of QCD matter. The same applies to the 
matter in the interior of neutron stars. The corresponding conservation of net 
baryon number is introduced into the grand canonical statistical model of 
QCD matter via the "baryo-chcmical potential" fiB- 

We employ this terminology to draw a phase diagram of QCD matter in 
figure 1, in the variables T and /zb- Note that /tts is high at low energies of 
collisions creating a matter fireball. In a head-on collision of two mass 200 
nuclei at -^i = 15GeV the fireball contains about equal numbers of newly 
created quark-antiquark pairs (of zero net baryon number), and of initial 
valence quarks. The accomodation of the latter, into created hadronic species, 
thus requires a formidable redistribution task of net baryon number, reflecting 
in a high value of /xb. Conversely, at LHC energy {^/s=5.5TeV in Pb-I-Pb 
collisions), the initial valence quarks constitute a mere 5% fraction of the total 
quark density, correspondingly requiring a small value oi hb- In the extreme, 
big bang matter evolves toward hadronization (at T=170 MeV) featuring a 
quark over antiquark density excess of 10~^ only, resulting in /i^ w 0. 

Note that the limits of existence of the hadronic phase are not only reached 
by temperature increase, to the so-called Hagedorn value Th (which coincides 
with Tcrit at ^B 0), but also by density increase to £i > (5 — 10) qq: "cold 
compression" beyond the nuclear matter ground state baryon density qq of 
about 0.16 Bj frr? . We are talking about the deep interior sections of neutron 
stars or about neutron star mergers [14, 15, 16]. A sketch of the present view 
of the QCD phase diagram [11, 12, 13] is given in Fig. 1. It is dominated by 
the parton-hadron phase transition line that interpolates smoothly between 
the extremes of predominant matter heating (high T, low /is) and predom- 
inant matter compression (T 0, \xb > 1 GeV). Onward from the latter 
conditions, the transition is expected to be of first order [17] until the critical 
point of QCD matter is reached at 200 < fiB (E) < 500 MeV. The relatively 
large position uncertainty reflects the preliminary character of Lattice QCD 
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calculations at finite ^.b [H, 12, 13]. Onward from the critical point, E, the 
phase transformation at lower is a cross-over [13]. 




(MeV) 

Fig. 1. Sketch of the QCD matter phase diagram in the plane of temperature T and 
baryo-chemical potential fj,B- The parton-hadron phase transition line from lattice 
QCD [10, 11, 12, 13] ends in a critical point E. A cross-over transition occurs at 
smaller ^lb ■ Also shown are the points of hadro-chemical freeze-out from the grand 
canonical statistical model. 



We note, however, that these estimates represent a major recent advance 
of QCD lattice theory which was, for two decades, believed to be restricted 
to the fj,B = situation. Onward from the critical point, toward lower /xs, 
the phase transformation should acquire the properties of a rapid cross-over 
[18], thus also including the case of primordial cosmological expansion. This 
would finally rule out former ideas, based on the picture of a violent first or- 
der " explosive" cosmological hadronization phase transition, that might have 
caused non-homogeneous conditions, prevailing during early nucleo-synthesis 
[19], and fluctuations of global matter distribution density that could have 
served as seedlings of galactic cluster formation [20]. However, it needs to be 
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stressed that the conjectured order of phase transformation, occuring along 
the parton - hadron phase boundary line, has not been unambiguously con- 
firmed by experiment, as of now. 

On the other hand, the position of the QCD phase boundary at low /zs 
has, in fact, been located by the hadronization points in the T, fis plane 
that are also illustrated in Fig. 1. They are obtained from statistical model 
analysis [21] of the various hadron multiplicities created in nucleus- nucleus 
collisions, which results in a [T, jis] determination at each incident energy, 
which ranges from SIS via AGS and SPS to RHIC energies, i.e. 3 < ^/s < 
200 GeV. Toward low /zb these hadronic freeze-out points merge with the 
lattice QCD parton-hadron coexistence line: hadron formation coincides with 
hadronic species freeze-out. These points also indicate the domain of the 
phase diagram which is accessible to relativistic nuclear collisions. The domain 
at /is > 1.5 GeV which is predicted to be in a further new phase of QCD 
featuring color-flavor locking and color superconductivity [22] will probably 
be accessible only to astrophysical observation. 

One may wonder how states and phases of matter in thcrmodynamical 
equilibrium - as implied by a description in grand canonical variables - can 
be sampled via the dynamical evolution of relativistic nuclear collisions. Em- 
ploying heavy nuclei, A w 200, as projectiles/targets or in colliding beams 
(RHIC, LHC), transverse dimensions of the primordial interaction volume do 
not exceed about 8/m, and strong interaction ceases after about 20 /m/c. We 
note, for now, that the time and dimension scale of primordial perturbative 
QCD interaction at the microscopic partonic level amounts to subfractions 
of 1 fm/c, the latter scale, however, being representative of non perturbative 
processes (confinement, "string" formation etc.). The A+A fireball size thus 
exceeds, by far, the elementary non perturbative scale. An equilibrium quark 
gluon plasma represents an extended non-perturbative QCD object, and the 
question whether its relaxation time scale can be provided by the expan- 
sion time scale of an A-l-A collision, needs careful examination. Reassuringly, 
however, the hadrons that are supposedly created from such a preceding non- 
perturbative QGP phase at top SPS and RHIC energy, do in fact exhibit 
perfect hydrodynamic and hadrochemical equilibrium, the derived [T, hb] 
values [21] thus legitimately appearing in the phase diagram, Fig. 1. 

In the present book we will order the physics obscrvablcs to be treated, 
with regard to their origin from successive stages that characterize the overall 
dynamical evolution of a relativistic nucleus-nucleus collision. In rough outline 
this evolution can be seen to proceed in three major steps. An initial period 
of matter compression and heating occurs in the course of interpenetration 
of the projectile and target baryon density distributions. Inelastic processes 
occuring at the microscopic level convert initial beam longitudinal energy to 
new internal and transverse degrees of freedom, by breaking up the initial 
baryon structure functions. Their partons thus acquire virtual mass, popu- 
lating transverse phase space in the course of inelastic perturbative QCD 
shower multiplication. This stage should be far from thermal equilibrium, 
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initially. However, in step two, inelastic interaction between the two arising 
parton fields (opposing each other in longitudinal phase space) should lead 
to a pile-up of partonic energy density centered at mid-rapidity (the longitu- 
dinal coordinate of the overall center of mass). Due to this mutual stopping 
down of the initial target and projectile parton fragmentation showers, and 
from the concurrent decrease of parton virtuality (with decreasing average 
square momentum transfer Q'^) there results a slowdown of the time scales 
governing the dynamical evolution. Equilibrium could be approached here, 
the system "lands" on the T, /z plane of Fig. 1, at temperatures of about 300 
and 200 MeV at top RHIC and top SPS energy, respectively. The third step, 
system expansion and decay, thus occurs from well above the QCD parton- 
hadron boundary line. Hadrons and hadronic resonances then form, which 
decouple swiftly from further inelastic transmutation so that their yield ratios 
become stationary (" frozen-out" ) . A final expansion period dilutes the system 
to a degree such that strong interaction ceases all together. 

It is important to note that the above description, in terms of succes- 
sive global stages of evolution, is only valid at very high energy, e.g. at 
and above top RHIC energy of ^/s = 200GeV. At this energy the target- 
projectile interpenetration time 2ii/7 = 0.12fm/c, and thus the interpenetra- 
tion phase is over when the supposed next phase (perturbative QCD shower 
formation at the partonic level by primordial, "hard" parton scattering) set- 
tles, at about 0.25fm/c. "Hard" observables (heavy flavour production, jets, 
high pt hadrons) all originate from this primordial interaction phase. On the 
other hand it is important to realize that at top SPS energy, ^/s = 17.3GeV, 
global interpenetration takes as long as 1.5fm/c, much longer than microscopic 
shower formation time. There is thus no global, distinguishable phase of hard 
QCD mechanisms: they are convoluted with the much longer interpenetration 
time. During that it is thus impossible to consider a global physics of the 
interaction volume, or any equilibrium. Thus we can think of the dynami- 
cal evolution in terms of global " states" of the system's dynamical evolution 
(such as local or global equilibrium) only after about 2-3fm/c, just before 
bulk hadronization sets in. Whereas at RHIC, and even more ideally so at the 
LHC, the total interaction volume is "synchronized" at times below 0.5fm/c, 
such that a hydrodynamic description becomes possible: we can expect that 
"flow" of partons sets in at this time, characterized by extremely high parton 
density. The dynamics at such early time can thus be accessed in well deflned 
variables (e.g. elliptic flow or jet quenching). 

In order to verify in detail this qualitative overall model, and to ascertain 
the existence (and to study the properties) of the different states of QCD that 
are populated in sequence, one seeks observable physics quantities that convey 
information imprinted during distinct stages of the dynamical evolution, and 
"freezing-out" without significant obliteration by subsequent stages. Ordered 
in sequence of their formation in the course of the dynamics, the most relevant 
such observables are briefly characterized below: 
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1. Suppression of J/??^ and Y production by Debye- screening in the QGP. 
These vector mesons result from primordial, pQCD production of cc and 
bb pairs that would hadronize unimpeded in elementary collisions but 
are broken up if immersed into a npQCD deconfined QGP, at certain 
characteristic temperature thresholds. 

2. Suppression of dijets which arise from primordial qq pair production frag- 
menting into partonic showers (jets) in vacuum but being attenuated by 
QGP-medium induced gluonic bremsstrahlung: Jet quenching in A-|-A 
collisions. 

a) A variant of this: any primordial hard parton suffers a high, specific 
loss of energy when traversing a deconfined medium: High px suppres- 
sion in A+A collisions. 

3. Hydrodynamic collective motion develops with the onset of (local) thermal 
equilibrium. It is created by partonic pressure gradients that reflect the 
initial collisional impact geometry via non-isotropies in particle emission 
called "directed" and "elliptic" flow. The latter reveals properties of the 
QGP, seen here as an ideal partonic fluid. 

a) Radial hydrodynamical expansion flow (" Hubble expansion" ) is a vari- 
ant of the above that occurs in central, head on collisions with cylinder 
symmetry, as a consequence of an iscntropic expansion. It should be 
sensitive to the mixed phase conditions characteristic of a first order 
parton-hadron phase transition. 

4. Hadronic "chemical" freeze-out fixes the abundance ratios of the hadronic 
species into an equilibrium distribution. Occuring very close to, or at 
hadronization, it reveals the dynamical evolution path in the [T, ^b] plane 
and determines the critical temperature and density of QCD. The yield 
distributions in A-|-A collisions show a dramatic strangeness enhancement 
eff'ect, characteristic of an extended QCD medium. 

5. Fluctuations, from one collision event to another (and even within a single 
given event) can be quantified in A-|-A collisions due to the high charged 
hadron multiplicity density (of up to 600 per rapidity unit at top RHIC 
energy). Such evcnt-by-evcnt (ebyc) fluctuations of pion rapidity density 
and mean transverse momentum (event "temperature"), as well as event- 
wise fluctuations of the strange to non-strange hadron abundance ratio 
(may) reflect the existence and position of the conjectured critical point 
of QCD (Fig. 1). 

6. Two particle Bose-Einstein-Correlations are the analog of the Hanbury- 
Brown, Twiss (HBT) effect of quantum optics. They result from the last 
interaction experienced by mesons, i.e. from the global decoupling stage. 
Owing to a near isentropic hadronic expansion they reveal information on 
the overall space-time-development of the "fireball" evolution. 

In an overall view the first group of obscrvablcs (1 to 2a) is anchored in 
established pQCD physics that is well known from theoretical and experimen- 
tal analysis of elementary collisions (e+e~ annihilation, pp and pp data). In 
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fact, the first generation of high baryon collisions, occuring at the micro- 
scopic level in A+A collisions, should closely resemble such processes. How- 
ever, their primary partonic products do not escape into pQCD vacuum but 
get attenuated by interaction with the concurrently developing extended high 
density medium , thus serving as diagnostic tracer probes of that state. The 
remaining observables capture snapshots of the bulk matter medium itself. 
After initial equilibration we may confront elliptic flow data with QCD dur- 
ing the corresponding partonic phase of the dynamical evolution employing 
thermodynamic [23] and hydrodynamic [24] models of a high temperature par- 
ton plasma. The hydro-model stays applicable well into the hadronic phase. 
Hadron formation (confinement) occurs in between these phases (at about 5 
microseconds time in the cosmological evolution). In fact relativistic nuclear 
collision data may help to finally pin down the mechanism(s) of this fascinat- 
ing QCD process [25, 26, 27] as we can vary the conditions of its occurence, 
along the parton-hadron phase separation line of Fig. 1, by proper choice of 
coUisional energy ^/s, and system size A, while maintaining the overall condi- 
tions of an extended imbedding medium of high energy density within which 
various patterns [11, 12, 13, 17, 18] of the hadronization phase transition may 
establish. The remaining physics observables (3a, 5 and 6 above) essentially 
provide for auxiliary information about the bulk matter system as it traverses 
(and emerges from) the hadronization stage, with special emphasis placed on 
manifestations of the conjectured critical point. 

The observables from 1 to 4 above will all be treated, in detail, in this 
book. We shall focus here on the bulk matter expansion processes of the 
primordially formed collisional volume, as reflected globally in the population 
patterns of transverse and longitudinal (rapidity) phase space (Section 3), and 
on the transition from partons to hadrons and on hadronic hadro- chemical 
decoupling, resulting in the observed abundance systematics of the hadronic 
species (Section 4). These Sections will be preceded by a detailed recapitulation 
of relativistic kinematics, notably rapidity, to which we shall turn now. 
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2 Relativistic Kinematics 

2.1 Description of Nucleus-Nucleus Collisions in terms of 
Light-Cone Variables 

In relativistic nucleus-nucleus collisions, it is convenient to use kinematic vari- 
ables which take simple forms under Lorentz transformations for the change 

of frame of reference. A few of them arc the light cone variables x+ and x-, 
the rapidity and pseudorapidity variables, y and tj. A particle is characterized 
by its 4-momentum, = p). In fixed target and collider experiments 
where the beam(s) define reference frames, boosted along their direction, it is 
important to express the 4-momentum in terms of more practical kinematic 
variables. 

Figure 2 shows the collision of two Lorentz contracted nuclei approaching 
each other with velocities nearly equal to velocity of light. The vertical axis 
represents the time direction with the lower half representing time before 
the collision and the upper half, time after the collision. The horizontal axis 
represents the spatial direction. Both the nuclei collide at (t, z) = (0, 0) and 
then the created fireball expands in time going through various processes 
till the created particles freeze-out and reach the detectors. The lines where 
t^ — = (note that \/t^ — = t, r being the proper time of the particle) 
along the path of the colliding nuclei define the light cone. The upper part of 
the light-cone, where — > 0, is the time-like region. In nucleus-nucleus 
collisions, particle production occurs in the upper half of the (t, z)-plane within 
the light-cone. The region outside the light cone for which — z"^ < is 
called space-like region. The space-time rapidity is defined as 



It could be seen that rjg is not defined in the space-like region. It takes the 
value of positive and negative infinity along the beam directions for which 
t = ±z respectively. A particle is "light- like" along the beam direction. Inside 
the light-cone which is time-like, rjs is properly defined. 

For a particle with 4-momentum p {po,Pt,Pz), the light-cone momenta 
are defined by 



p+ is called ''forward light-cone momentum" and p- is called "backward light- 
cone m,omentum^^ . 

For a particle traveling along the beam direction, has higher value of forward 
light-cone momentum and traveling opposite to the beam direction has lower 
value of forward light-cone momentum. The advantages of using light-cone 
variables to study particle production are the following. 




(1) 



P- =P0- Pz 



p+ =po+p^ 



(2) 
(3) 



10 



Michael Kliemant, Raghunath Sahoo, Tim Schuster, and Reinhard Stock 




Fig. 2. Description of heavy- ion collisions in one space (z) and one time (t) dimen- 
sion. 

1. The forward light-cone momentum of any particle in one frame is related 
to the forward light-cone momentum of the same particle in another boosted 
Lorentz frame by a constant factor. 

2. Hence, if a daughter particle c is fragmenting from a parent particle b, then 
the ratio of the forward light-cone momentum of c relative to that of b is 
independent of the Lorentz frame. 

Define 



Po+P'z 



pI 



P". 



(4) 



The forward light-cone variable x+ is always positive because c+ can't be 
greater than 6+. Hence the upper limit of a;+ is 1. x+ is Lorentz invariant. 
3. The Lorentz invariance of x+ provides a tool to measure the momentum of 
any particle in the scale of the momentum of any reference particle. 



2.2 The Rapidity Variable 

The co-ordinates along the beam line (conventionally along the z-axis) is 
called longitudinal and perpendicular to it is called transverse (x-y). The 3- 
momentum can be decomposed into the longitudinal (pz) and the transverse 
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(Pt), Pt being a vector quantity which is invariant under a Lorentz boost 
along the longitudinal direction. The variable rapidity "?/" is defined by 

y=oini — — . (5) 



2 \E-p, 

It is a dimensionless quantity related to the ratio of forward light-cone to 
backward light-cone momentum. The rapidity changes by an additive constant 
under longitudinal Lorentz boosts. 

For a free particle which is on the mass shell (for which E"^ = + w?), 
the 4-momentum has only three degrees of freedom and can be represented 
by (y, Pt)- {E,tpt) could be expressed in terms of (y, Pt) as 

E = rriT cosh y (6) 
Pz = Tot sink y (7) 

Tot being the transverse mass which is defined as 

TOy = TO^ -h Py. (8) 

The advantage of rapidity variable is that the shape of the rapidity distribution 
remains unchanged under a longitudinal Lorentz boost. When we go from 
CMS to LS, the rapidity distribution is the same, with the y-scale shifted by 
an amount equal to ycm- This is shown below. 



Rapidity of Center of Mass in the Laboratory System 



The total energy in the CMS system is Ecm = The energy and momentum 
of the CMS in the LS are jcmVs and PcmlcmVs respectively. The rapidity of 
the CMS in the LS is 



= — In 



— In 
2 



I -Per 



(9) 



It is a constant for a particular Lorentz transformation. 



Relationship between Rapidity of a particle in LS and rapidity in 
CMS 

The rapidities of a particle in the LS and CMS of the collision are respectively, 
y =\ln ^-§3^^ and y* = \ln { ^,^pl ^ ■ Inverse Lorentz transformations on 
E and Pz give 
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y=^ln 

= - In 
2 

y = y* + Vc 







IE" + , 




[^(E*+fipl)-^(pE*+pl) 


\E*+P*A 


+ lln 


'l + f3' 




[e*-p%\ 







(10) 

(11) 



Hence the rapidity of a particle in the laboratory system is equal to the sum 
of the rapidity of the particle in the center of mass system and the rapidity 
of the center of mass in the laboratory system. It can also be state that the 
rapidity of a particle in a moving (boosted) frame is equal to the rapidity in its 
own rest frame minus the rapidity of the moving frame. In the non-relativistic 
limit, this is like the subtraction of velocity of the moving frame. However, 
this is not surprising because, non-relativistically, the rapidity y is equal to 
longitudinal velocity /?. Rapidity is a relativistic measure of the velocity. This 
simple property of the rapidity variable under Lorentz transformation makes 
it a suitable choice to describe the dynamics of relativistic particles. 



Relationship between Rapidity and Velocity 



Consider a particle traveling in 2:-direction with a longitudinal velocity f3. The 
energy E and the longitudinal momentum p^ of the particle are 



E 

Pz 



7m 
7/3m 



(12) 
(13) 



where m is the rest mass of the particle. Hence the rapidity of the particle 
traveling in ^-direction with velocity /? is 



y/3 



— In 
2 

= - In 
2 



E+p, 



E-p, 
1 + /3' 
1-/3, 



\ln 



7m + 7/3m 



7m — 7/3m 



(14) 



Note here that yfj is independent of particle mass. In the non-relativistic limit 
when (3 is small, expanding in terms of /? leads to 

2/;3=/3+0(/33) (15) 
Thus the rapidity of the particle is the relativistic realization of its velocity. 



Beam Rapidity 

We know, 

E = m,T cosh y, pz ~ mx sink y and nt^ = m? + p|.. 
For the beam particles, = 0. 
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Hence, E = rrn, cosh yi, and = mi, sink yi,, 

where m{, and yi, are the rest mass and rapidity of the beam particles. 



yb = cosh ^ {E/rrib) 
\/snn 



= cosh ^ 
and 



2 m„ 



(16) 



yt = sink ^ {pz/rrih) (17) 
Here m„ is the mass of the nucleon. Note that the beam energy E = y'sjvAr/2. 

Example 1. For the nucleon- nucleon center of mass energy ^smn =9.1 GeV, 
the beam rapidity y^ = cosh~^ ( 2x0 oss ) ~ ^'^^ 
For p-hp collisions with lab momentum 100 GeV/c, 
y, = s^nh-^ (^) = sznh-^ {^) = 5.36 

and for Pb+Pb collisions at SPS with lab energy 158 AGeV, yb = 2.92. 

Rapidity of the CMS in terms of Projectile and Target Rapidities 

Let us consider the beam particle "6" and the target particle "o" . 

bz = rriT sink yb = rub sink yb- This is because pt of beam particles is 

zero. Hence 

yb = sinh"^ {h^/mb). (18) 

The energy of the beam particle in the laboratory frame is 
60 = rriT cosh yb = rub cosh yb- 

Assuming target particle a has longitudinal momentum Uz, its rapidity in the 
laboratory frame is given by 

ya = sinh~'^ {az/nia) (19) 

and its energy 

ao = ma cosh ya- (20) 

The CMS is obtained by boosting the LS by a velocity of the center-of-mass 
frame /3cm such that the longitudinal momentum of the beam particle 6* 
and of the target particle a* are equal and opposite. Hence /3cm satisfies the 
condition, 

a*z = lcm{az-f3cmao) = -b*z = - 7cm (&z- /3cm &o), where 7cm = yxT^- 
Hence, 

/3cm = (21) 

ao + 60 

We know the rapidity of the center of mass is 
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I In 



"l+/?c 
1-Pc 



Using equations 21 and 22, wc get 



Vc 



- In 
2 



ao 



ao- Uz + bo 



(22) 



(23) 



Writing energies and momenta in terms of rapidity variables in the LS, 
nia cosh Da + rua sink ya + mt cosh yb + mt sinh yb 



Vcm — 2 



rua cosh ya — rua sinh ya + rrib cosh yb — rub sinh yb 



1/ X 1 , 



nia e*'" + rub 



nia e*"' + nib e^" 
For a symmetric collision (for m-a = nib), 



= -^{ya + yb) 



Rapidities of a and b in the CMS are 



yl = ya-ycm = --^{yb-ya) 



Vb 



Vb - Vc 



liVb 



ya)- 



(24) 



(25) 



(26) 

(27) 



Given the incident energy, the rapidity of projectile particles and the rapidity 
of the target particles can thus be determined. The greater is the incident 
energy, the greater is the separation between the projectile and target rapidity. 
Central Rapidity The region of rapidity mid-way between the projectile and 
target rapidities is called central rapidity. 



Example 2. In p+p collisions at a laboratory momentum of 100 GeV/c, beam 
rapidity yb = 5.36, target rapidity j/a = and the central rapidity « 2.7. 

Mid-rapidity in Fixed target and Collider Experiments 

In fixed-target experiments (LS), ytarget = 0. 

yiab = ytarget + yprojectile = ybeam Hencc mid-rapidity in fixed-target 
experiment is given by, 

= ybeam/2. (28) 

In collider experiments (center of mass system), 

yprojectile — ytarget yCMS ybeam f'^- 

Hence, mid-rapidity in CMS system is given by 

ymf/ = {yprojectile + ytarget) l'^. = 0. (29) 

This is valid for a symmetric energy collider. The rapidity difference is given 
by yprojectile - ytarget = "^VCMS and this iucreascs with energy for a collider 
as y increases with energy. 
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Light-cone variables and Rapidity 

Consider a particle having rapidity y and the beam rapidity is yi,. The particle 
has forward light-cone variable x+ with respect to the beam particle 

^ _ Po+Pz 
^ P^o+pI 

= ^^e^-y" (30) 

where is the transverse mass of c. Note that the transverse momentum of 
the beam particle is zero. Hence, 

y = yb + /n a;+ + in [ ^ J (31) 

Similarly, relative to the target particle a with a target rapidity t/o, the back- 
ward light-cone variable of the detected particle c is a;_. X- is related to y 
by 

(32) 



TO. 



and conversely. 



TO, 



y=-Va-lnx_-ln[^\. (33) 



In general, the rapidity of a particle is related to its light-cone momenta by 

V=\l^ f— ) (34) 



2 \V 

Note that in situations where there is a frequent need to work with boosts 
along z-direction, it's better to use (y, px) for a particle rather than using it's 
3-momentum, because of the simple transformation rules for y and pt under 
Lorentz boosts. 



2.3 The Pseudorapidity Variable 



Let us assume that a particle is emitted at an angle Q relative to the beam 
axis. Then its rapidity can be written as 

y = hln(^) = \ ln 



2 ""' yE-PL j 
and hence 



^ m^+p^+p cos 9 
'\Jm'^+p^ —p cos 



. At very high energy, p ^ m 



- In 
2 



p + p cos 



P — P COS f 

In tan 0/2 = rj 



(35) 



r] is called the pseudorapidity. Hence at very high energy. 
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y Ki T] = —In tan 9/2. 



(36) 



In terms of the momentum, rj can be re-written as 



V =lln 



Pz 



(37) 



6 is the only quantity to be measured for the determination of pseudorapid- 
ity, independent of any particle identification mechanism. Pseudorapidity is 
defined for any value of mass, momentum and energy of the collision. This 
also could be measured with or without momentum information which needs 
a magnetic field. 

Change of variables from (y, p^) to (r;, pt) 



-i — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — p 
• 0-e% Au+Au200GeV 25 35% 




Fig. 3. The mid-rapidity dNch/dr} for Au+Au collisions at v'ijvjv = 200 GeV [28]. 



By equation 37, 



= J\p+I^ (38) 

\ m-Pz 



= (39) 

Adding both of the equations, we get 

|p| = "Pt cosh T] (40) 
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Pt = v|pp— pf . By subtracting the above equations, we get 

Pz = Pt sink rj 
Using these equations in the definition of rapidity, we get 
1 



y 



Ir 



\/Px cosh? 'q + w? + Pt sink rj 
\/p\ cosh? r] + m'^ — pr sinh rj 



Similarly r] could be expressed in terms of y as, 



sj m?p cosh^ y — m? + niT sinh y 
rriij. cosh? y — w? — mr sinh y 



(41) 



(42) 



(43) 



The distribution of particles as a function of rapidity is related to the distri- 
bution as a function of pseudorapidity by the formula 



dN 
dridpj 



dN 



cosh'^ y dydpT 



(44) 



In the region y ^ 0, the pseudorapidity distribution {dN/drj) and the rapidity 
distribution (dN/dy) which are essentially the pr-integrated values of 

and ^y^^ respectively, are approximately the same. In the region y ~ 0, there 
is a small ^'depression" in dN/dr} distribution compared to dN/dy distribu- 
tion due to the above transformation. At very high energies where dN/dy 
has a mid-rapidity plateau, this transformation gives a small dip in dN/drj 
around 77 « (see Figure 3). However, for a massless particle like photon, 
the dip in dN/dr] is not expected (which is clear from the above equation). 
Independent of the frame of reference where 77 is measured, the difference in 
the maximum magnitude of dN/ drj appears due to the above transformation. 
In the CMS, the maximum of the distribution is located at j/ « « and the 
77-distribution is suppressed by a factor -^/l — m? / < rn^ > with reference to 
the rapidity distribution. In the laboratory frame, however the maximum is 
located around half of the beam rapidity ry « t/f,/2 and the suppression factor 
is -^/l — rr? / < m^, > cosh? (j/5/2), which is about unity. Given the fact that 
the shape of the rapidity distribution is independent of frame of reference, the 
peak value of the pseudorapidity distribution in the CMS frame is lower than 
its value in LS. This suppression factor at SPS energies is 0.8 — 0.9. 



2.4 The Invariant Yield 

First we show ^ is Lorentz invariant. The differential of Lorentz boost in 
longitudinal direction is given by 



dpi = i{dp,-fidE). 



(45) 
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Taking the derivative of the equation E"^ = -\- mP, we get 

EdE = p^dp,. (46) 
Using equations 45 and 46 we get 

dpi = -f{dp^ - P^^Y^) 

= ^E*. (47) 
E 

As pt is Lorentz invariant, multiplying px on both the sides and re-arranging 

gives 

rfV ^ d'^PT dp^ ^ d^ 

E* E E ' ^ ' 

In terms of experimentally measurable quantities, ^ could be expressed as 

— = dpT dy 

= pTdpTd(f)dy (49) 
= mTdniTdcjydy. (50) 

The Lorentz invariant differential cross-section = is the invariant 

yield. In terms of experimentally measurable quantities this could be expressed 
as 

Ed^a 1 d^N 



dp^ rriT dmTd(j)dy 
1 d^N 



27r rriT dmrdy 
1 d'^N 



2n Pt dprdy 



(51) 



To measure the invariant yields of identified particles equation 51 is used 
experimentally. 



2.5 Inclusive Production of Particles and the Feynman Scaling 
variable xp 

In a reaction of type 

beam + target — > A + anything 

where A is known is called an ^Hnclusive reaction" . The cross-section for par- 
ticle production could be written separately as functions of px and pl- 



o- = f{PT)g{PL)- 



(52) 
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This factorization is empirical and convenient because each of these factors 
has simple parameterizations which fit well to experimental data. 

Similarly the differential cross-section could be expressed by 

(Pa _ (Pa da 
dp^ pI dpL 



Define the variable 



[max) 



2Pl 



/s 



(55) 



a;i? is called the Feynman scaling variable: longitudinal component of the cross- 
section when measured in CMS of the collision, would scale i.e. would not 
depend on the energy ^/s. Instead of is measured which wouldn't 

depend on energy of the reaction, ^/s. This Feynman's assumption is valid 
approximately. 

The differential cross-section for the inclusive production of a particle is 
then written as 

d£^^ = ('^^ 

Feynman's assumption that at high energies the function F{s, xp, Pt) be- 
comes asymptotically independent of the energy means: 
lims^ooF{s, xf, Pt) = F{xf, Pt) = /(pt) gixp) 



2.6 The PT-Distribution 

The distribution of particles as a function of pr is called pT-distribution. 
Mathematically, 

^ = (57) 
dpT 2tt IptMIptI 

whore dN is the number of particles in a particular pj^-bin. People usually 
plot p^l^p^ as a function of px taking out the factor l/27r which is a con- 
stant. Here px is a scalar quantity. The low-p^ part of the pi -spcctrum is 
well described by an exponential function having thermal origin. However, to 
describe the whole range of the pT, one uses the Levy function which has an 
exponential part to describe low-pT and a power-law function to describe the 
hight-pT part which is dominated by hard scatterings (high momentum trans- 
fer at early times of the collision). The inverse slope parameter of pT-spectra 
is called the effective temperature (Tg//), which has a thermal contribution 
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because of the random kinetic motion of the produced particles and a con- 
tribution from the collective motion of the particles. This will be described 
in details in the section of freeze-out properties and how to determine the 
chemical and kinetic freeze-out temperatures experimentally. 

The most important parameter is then the mean pT which carries the in- 
formation of the effective temperature of the system. Experimentally, {pr) is 
studied as a function of which is the measure of the entropy density of 
the system. This is like studying the temperature as a function of entropy 
to see the signal of phase transition. The phase transition is of 1st order if a 
plateau is observed in the spectrum signaling the existence of latent heat of 
the system. This was first proposed by L. Van Hove [29]. 

The average of any quantity A following a particular probability distribu- 
tion f{A) can be written as 

_ SAmdA 

- Jf{A)dA ■ ^^^^ 



Similarly, 



ij^) dpr 

Cpt dpT Pri^i^) 
J^PT dpTij^) 

PT dpT PT Upt) 
PT dpT,f{pT) 



(59) 



where 27r pt dpT is the phase space factor and the pr-distribution function is 
given by 

dN dN 
fipr) = -1 — = — (60) 

dpT pTdpT 

Example 3. Experimental data on py-spectra are sometimes fitted to the ex- 
ponential Boltzmann type function given by 

f{p^) = J-^ ^ Ce-™-/^»//. (61) 
Pt dpT 

The (mr) could be obtained by 

Jf^PT dpT rriT exp.{-mT/Teff) 



{ttit) = 



Jo Pt dpT exp.{-mT/Teff) 
2T^fj: + 2maT^ff +ml 
mo + Teff 



(62) 
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where mo is the rest mass of the particle. It can be seen from the above 
expression that for a massless particle 

{mr) = (pt) = "iTeff. (63) 

This also satisfies the principle of equipartition of energy which is expected 
for a massless Boltzmann gas in equilibrium. 

However, in experiments the higher limit of pt is a finite quantity. In that 
case the integration will involve an incomplete gamma function. 

2.7 Energy in CMS and LS 

For Symmetric Collisions {A + A) 

Consider the collision of two particles. In LS, the projectile with momentum 
Pi, energy Ei and mass rrii collides with a particle of mass m2 at rest. The 
4-momenta of the particles are 
Pi = (-E^i.Pi), P2 = (m2,0) 

In CMS, the momenta of both the particles are equal and opposite, the 4- 

momenta are 

pl = iEl,pl), pl = [El-vD 

The total 4-momentum of the system is a conserved quantity in the collision. 
In CMS, 

{.PI+P2? = {E^+E2f-{pi+P2? 

= {E^+E^f = El^ = 8. 

a/s is the total energy in the CMS which is the invariant mass of the CMS. 
In LS, 

(pi+P2)^ = ml + ml + 2Eim2. 
Hence 

Ecm = \fs = ^Jm\+m% + 2Eprojm2 (64) 

where Ei = Eproj, the projectile energy in LS. Hence it is evident here that 
the CM frame with an invariant mass y/s moves in the laboratory in the 
direction of pi with a velocity corresponding to: 
Lorentz factor, 

7cm = 1= — (65) 

^^s=^, (66) 

this is because E = 7m and 

Vcm = COSh~'^ 'Jem- (67) 
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Note 1. We know 

s = El^ = ml + ml + 2{Ei+E2 + Pi.p2) (68) 

For a head-on collision with TOi, m2 ^ Ei,E2 

- "^E.E^ (69) 

For two beams crossing at an angle 9, 

El^ = 2EiE2{l + cos 6) (70) 

The CM energy available in a collider with equal energies (E) for new particle 
production rises linearly with E i.e. 

Ecm ^ 2E (71) 

For a fixed-target experiment the CM energy rises as the square root of the 
incident energy: 

Ecm ^ V2m2Ei (72) 

Hence the highest energy available for new particle production is achieved at 
collider experiments. For example, at SPS fixed-target experiment to achieve 
a CM energy of 17.3 AGeV the required incident beam energy is 158 AGeV. 

Note 2. Most of the times the energy of the collision is expressed in terms of 
nucleon-nuclcon center of mass energy. In the nucleon-nucleon CM frame, two 
nuclei approach each other with the same boost factor 7. The nucleon-nucleon 
CM is denoted by ^/SNN and is related to the total CM energy by 

^ = A (73) 

This is for a symmetric collision with number of nucleons in each nuclei as 
A. The colliding nucleons approach each other with energy ^/snn/'^ and with 
equal and opposite momenta. The rapidity of the nucleon-nucleon center of 
mass is 

Vnn — and taking mi = m2 = rriN, the projectile and target nucleons are 
at equal and opposite rapidities. 

— _ — V^^^ _ (74^ 

Vproj — Vtarget — COSfl — VbeaTn' v'^/ 

Note 3. : Lorentz Factor 



M 2A ruN 

A y'SjVJV _ ^/SNN 

2A niN 2 mpf 

T?CMS 

beam (7^\ 

rriN 
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where E and M are Energy and Mass in CMS respectively. Assuming mass of 
the nucleon m^v ~ 1 GeV, the Lorentz factor is of the order of beam energy 
in CMS for a symmetric colhsion. 



For Asymmetric Collisions {A + B) 

During the early phase of relativistic nuclear collision research, the projectile 
mass was limited by accelerator-technical conditions (^*Ar at the Bcvalac, ^^Si 
at the AGS, ^^S at the SPS). Nevertheless, collisions with mass « 200 nuclear 
targets were investigated. Analysis of such collisions is faced with the problem 
of determining an " effective" center of mass frame, to be evaluated from the 
numbers of projectile and target participant nucleons, respectively. Their ratio 
- an thus the effective CM rapidity - depends on impact parameter. Moreover, 
this effective CM frame refers to soft hadron production only, whereas hard 
processes are still referred to the frame of nucleon-nucleon collisions. The light 
projectile on heavy target kinematics are described in [30] . 



2.8 Luminosity 

The luminosity is an important parameter in collision experiments. The reac- 
tion rate in a collider is given by 



R = aL (76) 



where, 

a = interaction cross-section 
L = luminosity (in cm^'^s^^) 



fn ^ (77) 



where, 

/ = revolution frequency 

A'^i, N2 = number of particles in each bunch 

n = number of bunches in one beam in the storage ring 

A = cross-sectional area of the beams 

L is larger if the beams have small cross-sectional area. 



2.9 Collision Centrality 

In a collision of two nuclei, the impact parameter (&) can carry values from 
to i?i + i?2, where Ri and R2 are the diameters of the two nuclei. When 
6 = 0, it is called head-on collision. When collisions with < 6 < {Ri + R2) 
are allowed, it is called minimum-bias collision. In heavy-ion collisions, initial 
geometric quantities such as impact parameter and the collision geometry can 
not be directly measured experimentally. Contrary, it is however possible to 



24 



Michael Kliemant, Raghunath Sahoo, Tim Schuster, and Reinhard Stock 



12 10 & 6 i 2 <b ifm\> 

1 F~l 1 1 1 1 1 1 




Fig. 4. A cartoon showing the centrality definition from the final state particle 
multiplicity and its correlation with the impact parameter (fe) and the number of 
participating nucleons (Npart) in the collisions. 

relate the particle multiplicity, transverse energy and the number of spectac- 
tor nucleons (measured by a "zero-degree calorimeter" ZDC) to the centrality 
of the collisions. 

It is straight forward to assume that on the average, 

1) the energy released in a collision is proportional to the number of nucleons 
participating in the collisions. 

2) the particle multiplicity is proportional to the participating nuclcon num- 
ber. 

Hence the particle multiplicity is proportional to the energy released in the col- 
lision. One can measure the particle multiplicity distribution or the transverse 
energy (Ex) distribution for minimum-bias collisions. Here the high values of 
particle multiplicity or Et correspond to central collisions and lower values 
correspond to more peripheral collisions. Hence the minimum-bias Ej- or mul- 
tiplicity distribution could be used for centrality determination in an collision 
experiment. Figure 4 shows the minimum-bias multiplicity (Nch) distribution 
used for the selection of collision centrality. The minimum-bias yield has been 
cut into successive intervals starting from the maximum value of Nch- The 
first 5% of the high Nch events correspond to top 5% central collisions. The 
correlation of centrality and the impact parameter with the number of partic- 
ipating nucleons has also been elaborated, in detail, by Glauber-type Monte 
Carlo calculations employing Woods-Saxon nuclear density distributions. 
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2.10 Number of Participants and Number of Binary Collisions 

Experimentally there is no direct way to estimate the number of participating 
nucleons {Npart) and the number of binary collisions (A^bm) in any event, for 
an given impact parameter. The Glauber model calculation is performed to 
estimate the above two quantities as a function of the impact parameter. The 
Glauber model treats a nucleus-nucleus collision as a superposition of many 
independent nucleon-nucleon (N — N) collisions. This model depends on the 
nuclear density profile (Woods-Saxon) and the non-diffractive inelastic N + N 
cross-sections. The Woods-Saxon distribution is given by 



where, r is the radial distance from the center of the nucleus, tq is the mean 
radius of the nucleus, c is the skin depth of the nucleus and po is the nuclear 
density constant. The parameters tq and c are measured in electron- nucleus 
scattering experiments, po is determined from the overall normalization con- 
dition 



where A is the mass number of the nucleus. 

There are two separate implementations of Glauber approach: Optical and 
Monte Carlo (MC). In the Optical Glauber approach, Npart and Ni,in are es- 
timated by an analytic integration of overlapping Woods-Saxon distributions. 

The MC Glauber calculation proceeds in two steps. First the nucleon po- 
sition in each nucleus is determined stochastically. Then the two nuclei are 
"collided", assuming the nucleons travel in a straight line along the beam 
axis (this is called eikonal approximation). The position of each nucleon in 
the nucleus is determined according to a probability density function which is 
typically taken to be uniform in azimuth and polar angles. The radial prob- 
ability function is modeled from the nuclear charge densities extracted from 
electron scattering experiments. A minimum inter-nucleon separation is as- 
sumed between the positions of nucleons in a nucleus, which is the charac- 
teristic length of the repulsive nucleon-nucleon force. Two colliding nuclei 
are simulated by distributing A nucleons of nucleus A and B nucleons of nu- 
cleus B in 3-dimensional co-ordinate system according to their nuclear density 
distribution. A random impact parameter b is chosen from the distribution 
dcr/db = 27r6. A nucleus- nucleus collision is treated as a sequence of indepen- 
dent nucleon-nucleon collisions with a collision taking place if their distance 
D in the transverse plane satisfies 



where o"^^ is the total inelastic nucleon-nucleon cross-section. An arbitrary 
number of such nucleus- nucleus collisions are performed by the monte carlo 




(79) 




(80) 
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and the resulting distributions of da/Npart and da/Nhin, da/db are deter- 
mined. Here Npart is defined as the total number of nuclcons that underwent 
at least one interaction and A^(,j„ is the total number of interactions in an 
event. These histograms are binned according to fractions of the total cross- 
sections. This determines the mean values oi Npart and Nf,in for each centrality 
class. The systematic uncertainties in these values are estimated by varying 
the Wood-Saxon parameters, by varying the value of cr^^i and from the un- 
certainty in the determination of total nucleus-nucleus cross-section. These 
sources of uncertainties are treated as fully correlated in the final systematic 
uncertainty in the above measured variables. 

When certain cross-sections scale with number of participants, those are 
said to be associated with "soft" processes: small momentum transfer pro- 
cesses. The low-pT hadron production which accounts for almost 95% of 
the bulk hadron multiplicity comes in the "soft processes" . These soft pro- 
cesses are described by phenomenological non-perturbative models. Whereas, 
in "hard" QCD processes like jets, charmonia, other heavy flavor and pro- 
cesses associated with high-pr phenomena, the cross-section scales with the 
number of primordial target/projectile parton collisions. This is estimated in 
the above Glauber formalism as the total number of inelastic participant- 
participant collisions. For the hard processes the interaction is at partonic 
level with large momentum transfer and is governed by pQCD. Ncoii is always 
higher than Npart- when Npart grows like A, NcoU grows like A^/^. 

Sometimes, to study the contribution of soft and hard processes to any 
cross-section one takes a two-component model like: 

cross - section = (1 - f) Npart + / NcoU (81) 
where / is the fractional contribution from hard processes. 
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3 Bulk Hadron Production in A+A Collisions 

We will now take an overall look at bulk hadron production in nucleus-nucleus 
collisions. In view of the high total cm. energies involved at e.g. top SPS 
{El'^ w 3.3 TeV) and top RHIC (38 TeV) energies, in central Pb+Pb (SPS) 
and Au-|-Au (RHIC) collisions, one can expect an extraordinarily high spatial 
density of produced particles. The average number of produced particles at 
SPS energies is w 1600, while at RHIC multiplicities of w 4000 arc reached. 
Thus, as an overall idea of analysis, one will try to relate the observed flow 
of energy into transverse and longitudinal phase space and particle species to 
the high energy density contained in the primordial interaction volume, thus 
to infer about its contained matter. 

Most of the particles under investigation correspond to "thermal" pions 
{pt up to 2 GeV) and, in general, such thermal hadrons make up for about 
95% of the observed multiplicity: the bulk of hadron production. Their dis- 
tributions in phase space will be illustrated in the subsections below. This 
will lead to a first insight into the overall reaction dynamics, and also set 
the stage for consideration of the rare signals, imbedded in this thermal bulk 
production: direct photons, jets, heavy flavors, which are the subject of later 
chapters in this volume. 

3.1 Particle Multiplicity and Transverse Energy Density 

Particle production can be assessed globally by the total created transverse 
energy, the overall result of the collisional creation of transverse momentum px 
or transverse mass (tot = ^/Pt + '^o)' microscopic Icivcl. Fig. 5 shows 

the distribution of total transverse energy Et = -^(^i) ' sinO resulting from 

i 

a calorimetric measurement of energy flow into calorimeter cells centered at 
angle 6i relative to the beam [31], for ^^S'-l-^^^Au collisions at ^/s = 20 GeV, 
and for 208p5+208pt, coUisions at = 17.3 GeV. 

The shape is characteristic of the impact parameter probability distribution 
(for equal size spheres in the Pb-|-Pb case) . The turnoff at Et = 520 GeV 
indicates the point where geometry runs out of steam, i.e. where 6 ^ 0, a con- 
figuration generally referred to as a "central collision". The adjacent shoulder 
resiilts from genuine event by event fluctuations of the actual number of par- 
ticipant nucleons from target and projectile (recall the diffuse Woods-Saxon 
nuclear density profiles), and from experimental factors like calorimeter reso- 
lution and limited acceptance. The latter covers 1.3 units of pseudo-rapidity 
and contains mid-rapidity r]mid = 2.9. Re-normalizing [31] to Zi 77 = 1 leads 
to dET / dri{mid) = AQQ GeV , in agreement with the corresponding WA80 
result [32]. Also, the total transverse energy of central Pb+Pb collisions at 
= 17.3 GeV turns out to be about 1.2 TeV. As the definition of a central 
collision, indicated in Fig. 5, can be shown [33] to correspond to an average 
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Fig. 5. Minimum bias distribution of total transverse energy in Pb+Pb collisions 
at ^/s = 17.3 GeV , and S+Au collisions at \fs = 20 GeV , in the rapidity interval 
2.1<y< 3.4, from [31]. 



nucleon participant number of Npart = 370 one finds an average total trans- 
verse energy per nucleon pair, of Et/ (0.5 Npart) = 6.5 GeV. After proper 
consideration of the baryon pair rest mass (not contained in the calorimet- 
ric Et response but in the corresponding yi) one concludes [31] that the 
observed total Et corresponds to about 0.6 i?™"^, the maximal Et derived 
from a situation of "complete stopping" in which the incident y/s gets fully 
transformed into internal excitation of a single, ideal isotropic fireball located 
at mid-rapidity. The remaining fraction of EJp"-^ thus stays in longitudinal 
motion, reflecting the onset, at SPS energy, of a transition from a central fire- 
ball to a longitudinally extended "firc-tubc", i.e. a cylindrical volume of high 
primordial energy density. In the limit of much higher ^ one may extrap- 
olate to the idealization of a boost invariant primordial interaction volume, 
introduced by Bjorkcn [34]. 

We shall show below (section 3.2) that the charged particle rapidity distri- 
butions, from top SPS to top RHIC energies, do in fact substantiate a devel- 
opment toward a boost-invariant situation. One may thus employ the Bjorkcn 
model for an estimate of the primordial spatial energy density e, related to 
the energy density in rapidity space via the relation [34] 

11 dET 
TTit^ To ay 

where the initially produced collision volume is considered as a cylinder of 

length dz — Tody and transverse radius R oc A^^^. Inserting for ttR"^ the lon- 
gitudinally projected overlap area of Pb nuclei colliding near head-on (" cen- 
trally"), and assuming that the evolution of primordial pQCD shower mul- 
tiplication (i.e. the energy transformation into internal degrees of freedom) 
proceeds at a time scale tq < 1/m/c, the above average transverse energy 
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density, of dEx/dy = 400 GeV at top SPS energy [31, 32] leads to the esti- 
mate 

e(To = 1/m) = 3.0 ± 0.6 GeV/fm^, (83) 

thus exceeding, by far, the estimate of the critical energy density cq obtained 
from lattice QCD (see below), of about 1.0 GeV/fm^. Increasing the collision 
energy to \/i = 2Q0GeV for Au+Au at RHIC, and keeping the same formation 
time, To = 1 fm/c (a conservative estimate as we shall show in section 3.4), 
the Bjorken estimate grows to e w 6.0 ± 1 GeV/ fm^. This statement is based 
on the increase of charged particle multiplicity density at mid-rapidity with 
■^i, as illustrated in Fig. 6. From top SPS to top RHIC energy [35] the density 
per participant nucleon pair almost doubles. However, at \/s = 200 GeV the 
formation or thermalization time tq, employed in the Bjorken model [34], 
was argued [36] to be shorter by a factor of about 4. We will return to such 
estimates of tq in section 3.5 but note, for now, that the above choice of 
To = 1 fm/c represents a conservative upper limit at RHIC energy. 
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Fig. 6. Charged liadron rapidity density at mid-rapidity vs. ^/s, compiled from 

e'^e^, pp, pp and A+A collisions [37]. 



These Bjorken-estimates of spatial transverse energy density are confronted 
in Fig. 7 with lattice QCD results obtained for three dynamical light quark 
flavors [38], and for zero baryo- chemical potential (as is realistic for RHIC 
energy and beyond but still remains a fair approximation at top SPS energy 
where /xs ~ 250 MeV). The energy density of an ideal, relativistic parton gas 
scales with the fourth power of the temperature. 



(84) 
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where g is related to the number of degrees of freedom. For an ideal gluon 
gas, g = 16 7r^/30; in an interacting system the effective g is smaller. The 
results of Fig. 7 show, in fact, that the Stefan-Boltzmann limit ess is not 
reached, due to non perturbative effects, even at four times the critical tem- 
perature Tc = f 70 MeV. The density e/T^ = 3 is seen to ascend steeply, 
within the interval Tc ± 25 MeV. At Tc the critical QCD energy density 
e = 0.6 - 1.0 GeV/fm^. Relating the thermal energy density with the Bjorken 
estimates discussed above, one arrives at an estimate of the initial temper- 
atures reached in nucleus-nucleus collisions, thus implying thermal partonic 
equilibrium to be accomplished at time scale tq (see section 3.5). For the SPS, 
RHIC and LHC energy domains this gives an initial temperature in the range 
190 < T^-P^ < 220 MeV, 220 < T^^^<^ < 400 MeV (assuming [36] that tq 
decreases to about 0.3 fm/c here) and T'"^^' > 600 MeV, respectively. From 
such estimates one tends to conclude that the immediate vicinity of the phase 
transformation is sampled at SPS energy, whereas the dynamical evolution at 
RHIC and LHC energies dives deeply into the " quark-gluon-plasma" domain 
of QCD. We shall return to a more critical discussion of such ascertations in 
section 3.5. 

One further aspect of the mid-rapidity charged particle densities per par- 
ticipant pair requires attention: the comparison with data from elementary 
collisions. Fig. 6 shows a compilation of pp, pp and e^e~ data covering the 
range from ISR to LEP and Tevatron energies. 




Fig. 7. Lattice QCD results at zero baryon potential for energy density e/T* versus 
T/Tc with three light quark flavors, compared to the Stefan-Boltzmann- hmit ess of 
an ideal quark-gluon gas [38]. 



Global Properties of Nucleus-Nucleus Collisions 



31 




PHOBOS Preliminary 



WHii+H 



200 GeV 



130GeV 



o 



^l^^iiAAi 19.6 GeV 



Z 10 



e e' 







o pp/pp 
• PHOBOS 







200 



400 



N 



part 



Fig. 8. The total number of charged hadrons per participant pair shown as a func- 
tion of Npart in Au+Au collisions at three RHIC energies [37]. 



The data from e+e~ represent dNch/dy, the rapidity density along the event 
thrust axis, calculated assuming the pion mass [39] (the difference between 
dN/ dy and dN/ drj can be ignored here) . Remarkably, they superimpose with 
the central A+A collision data, whereas pp and pp show similar slope but 
amount to only about 60% of the AA and e+e~ values. This difference be- 
tween e~^e~ annihilation to hadrons, and pp or pp hadro-production has been 
ascribed [40] to the characteristic leading particle effect of minimum bias 
hadron-hadron collisions which is absent in e^e~. It thus appears to be re- 
duced in AA collisions due to subsequent interaction of the leading parton 
with the oncoming thickness of the remaining target /projectile density distri- 
bution. This naturally leads to the scaling of total particle production with 
Npart that is illustrated in Fig. 8, for three RHIC energies and minimum bias 
Au+Au collisions; the close agreement with e+e^ annihilation data is obvious 
again. One might conclude that, analogously, the participating nucleons get 
"annihilated" at high -y/s, their net quantum number content being spread 
out over phase space (as we shall show in the next section). 

3.2 Rapidity Distributions 

Particle production number in A-t-A collisions depends globally on ^/s and 
collision centrality, and differentially on pr and rapidity y, for each particle 
species i. Integrating over pT results in the rapidity distribution dNi/dy. Parti- 
cle rapidity y = sinh~^pL/MT (where Mr = \/ m? -f pf^), requires mass iden- 
tification. If that is unknown one employs pseudo-rapidity (77 — — ln[tan(^/2)]) 
instead. This is also chosen if the joint rapidity distribution of several unre- 
solved particle species is considered: notably the charged hadron distribution. 
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dN /dri W=130GeV 




Fig. 9. Left panel: charged particle pseudo-rapidity distribution in pp collisions at 
^ = 540 GeF [41]. Right panel: same in RHIC Au+Au collisions at = 130 GeV 
at different centralities [42] . Closed lines represent fits with the color glass condensate 
model [43]. 



We show two examples in Fig. 9. The left panel illustrates charged particle 
production inpp collisions studied by UAl at y/s — biOGeV [41]. Whereas the 
minimum bias distribution (dots) exhibits the required symmetry about the 
center of mass coordinate, 77 = 0, the rapidity distribution corresponding to 
events in which a W boson was produced (histogram) features, both, a higher 
average charged particle yield, and an asymmetric shape. The former effect 
can be seen to reflect the expectation that the W production rate increases 
with the "centrality" of pp collisions, involving more primordial partons as 
the collisional overlap of the partonic density profiles gets larger, thus also 
increasing the overall, softer hadro-production rate. The asymmetry should 
result from a detector bias favoring W identification at negative rapidity: the 
transverse W energy, of about 100 Ge^ would locally deplete the energy store 
available for associated soft production. If correct, this interpretation suggests 
that the wide rapidity gap between target and projectile, arising at such high 
•^s, of width Ay 2ln {2^cm), makes it possible to define local sub-intervals 
of rapidity within which the species composition of produced particles varies. 

The right panel of Fig. 9 shows charged particle pseudo-rapidity density 
distributions for Au+Au collisions at y/s = 130 GeV measured by RHIC ex- 
periment PHOBOS [42] at three different collision centralities, from "central" 
(the 6% highest charged particle multiplicity events) to semi-peripheral (the 
corresponding 35-45% cut). We will turn to centrality selection in more detail 
below. Let us first remark that the slight dip at mid-rapidity and, moreover, 
the distribution shape in general, are common to pp and Au+Au. This is 
also the case for e^e~ annihilation as is shown in Fig. 10 which compares 
the ALEPH rapidity distribution along the mean px ("thrust") axis of jet 
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production in e+e" at ^/s = 200 GeV [39] with the scaled PHOBOS-RHIC 
distribution of central Au+Au at the same y'i [37]. Note that the mid-rapidity 
values contained in Figs. 9 and 10 have been employed already in Fig. 6, which 
showed the overall ^/s dependence of mid-rapidity charged particle produc- 
tion. What we concluded there was a perfect scaling of A+A with e+e" data 
at ^/s > 20 GeV and a 40% suppression of the corresponding pp, pp yields. 
We see here that this observation holds, semi-quantitatively, for the entire 
rapidity distributions. These are not ideally boost invariant at the energies 
considered here but one sees in dNch/dy a relatively smooth "plateau" region 
extending over | y \< 1.5 — 2.5. 

The production spectrum of charged hadrons is, by far, dominated by soft 
pions {pt < 1 GeV/c) which contribute about 85% of the total yield, both in 
elementary and nuclear collisions. The evolution of the tt^ rapidity distribu- 
tion with ^/s is illustrated in Fig. 11 for central Au-|-Au and Pb-|-Pb collisions 
from AGS via SPS to RHIC energy, 2.7 < ^/i < 200 GeV [44]. 




Fig. 10. Pseudo-rapidity distribution of charged liadrous produced in central 
Au-|-Au collisions at ^/s = 200 GeV compared with e'^e" data at similax energy. 
The former data normalized by Npartfi. Prom ref. [37]. 



At lower ^/s the distributions are well described by single Gaussian fits [44] 
with a{y) nearly linearly proportional to the total rapidity gap Ay oc ln^/s 
as shown in the right hand panel of Fig. 11. Also illustrated is the prediction 
of the schematic hydrodynamical model proposed by Landau [45], 



r^cclni^J (85) 
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which pictures hadron production in high ^/s pp cohisions to proceed via a 
dynamics of initial complete "stopping down" of the reactants matter/energy 
content in a mid-rapidity fireball that would then expand via 1-dimensional 
ideal hydrodynamics. Remarkably, this model that has always been considered 
a wildly extremal proposal falls rather close to the lower ^/s data for central 
A+A collisions but, as longitudinal phase space widens approaching boost in- 
variance we expect that the (non-Gaussian) width of the rapidity distribution 
grows linearly with the rapidity gap Ay. LHC data will finally confirm this 
expectation, but Figs. 9 to 11 clearly show the advent of boost invariance, 
aheady at Vs = 200 GeV. 

A short didactic aside: At low y/s the total rapidity gap Ay — 2 — 3 does 
closely resemble the total rapidity width obtained for a thermal pion velocity 
distribution at temperature T = 120 — 150 MeV, of a single mid-rapidity 
fireball, the y-distribution of which represents the longitudinal component 
according to the relation [21] 

dN , 2 2toT2 2r2 

— — oc (m T H — H —5—) exp \—m ■ coshy/T\ (86) 

dy coshy cosh'^y 

where m is the pion mass. Any model of preferentially longitudinal expansion 
of the pion emitting source, away from a trivial single central "completely 
stopped" fireball, can be significantly tested only once Ay > 3 which occurs 
upward from SPS energy. The agreement of the Landau model prediction with 
the data in Fig. 11 is thus fortuitous, below y/s w 10 GeV, as any created 
fireball occupies the entire rapidity gap with pions. 
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Fig. 11. Left panel: negative pion rapidity distributions in central Au+Au and 
Pb-I-Pb coUisions from AGS via SPS to RHIC energies [44]. Right panel: the Gaus- 
sian rapidity width of pions versus -^s, confronted by Landau model predictions 
(soUd line) [44]. 



The Landau model offers an extreme view of the mechanism of " stopping" , by 
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which the initial longitudinal energy of the projectile partons or nucleons is 
inelastically transferred to produced particles and redistributed in transverse 
and longitudinal phase space, of which we saw the total transverse fraction in 
Fig. 5. Obviously e+e~ annihilation to hadrons represents the extreme stop- 
ping situation. Hadronic and nuclear collisions offer the possibility to analyze 
the final distribution in phase space of their non-zero net quantum numbers, 
notably net baryon number. Figure 12 shows the net-proton rapidity dis- 
tribution (i.e. the proton rapidity distribution subtracted by the antiproton 
distribution) for central Pb+Pb/Au+Au collisions at AGS {^/s = 5.5 GeV), 
SPS < 17.3 GeV) and RHIC (Vs = 200 GeV) [46]. With increasing en- 
ergy wc sec a central (but non-Gaussian) peak developing into a double-hump 
structure that widens toward RHIC leaving a plateau about mid-rapidity. The 
RHIC-BRAHMS experiment acceptance for p, p identification does unfortu- 
nately not reach up to the beam fragmentation domain at jjp = 5.4 (nor does 
any other RHIC experiment) but only to y ~ 3.2, with the consequence that 
the major fraction of p"^* is not accounted for. However the mid-rapidity re- 
gion is by no means net baryon free. At SPS energy the NA49 acceptance 
covers the major part of the total rapidity gap, and we observe in detail a 
net p distribution shifted down from pp = 2.9 by an average rapidity shift 
[46] of (Sy) = 1.7. From Fig. 12 we infer that (Sy) can not scale linearly with 
t/p « ln{2jCM) ~ luy/s for ever - as it does up to top SPS energy where 
(Sy) = 0.58 yp [46]. Because extrapolating this relation to ^/s = 200 GeV 
would result in (Sy) =3.1, and with yp w 5.4 at this energy we would expect 
to observe a major fraction of net proton yield in the vicinity of y = 2.3 which 
is not the case. A saturation must thus occur in the (Sy) vs. \/i dependence. 

The re-distribution of net baryon density over longitudinal phase space is, of 
course, only partially captured by the net proton yield but a recent study [47] 
has shown that proper inclusion of neutron"^ and hyperon production data at 
SPS and RHIC energy scales up, of course, the dN/dy distributions of Fig. 12 
but leaves the peculiarities of their shapes essentially unchanged. As the net 
baryon rapidity density distribution should resemble the final valence quark 
distribution the Landau model is ruled out as the valence quarks are seen to 
be streaming from their initial position at beam rapidity toward mid-rapidity 
(not vice versa). It is remarkable, however, to see that some fraction gets 
transported very far, during the primordial partonic non-equilibrium phase. 
We shall turn to its theoretical description in section 3.4 but note, for now, 
that pp collisions studied at the CERN ISR [48] lead to a qualitatively similar 
net baryon rapidity distribution, albeit characterized by a smaller (Sy). 



^ Neutrons are not directly measured in the SPS and RHIC experiments but their 
production rate, relative to protons, reflects in the ratio of tritium to "^He pro- 
duction measured by NA49 [47], applying the isospin mirror symmetry of the 
corresponding nuclear wave functions. 
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Fig. 12. Net proton rapidity distributions in central Au+Au/Pb+Pb collisions at 
AGS, SPS and RHIC energies [46, 47]. 



The data described above suggest that the stopping mechanism universally 
resides in the primordial, first generation of collisions at the microscopic level. 
The rapidity distributions of charged particle multiplicity, transverse energy 
and valence quarks exhibit qualitatively similar shapes (which also evolve 
similarly with y^) in pp, pp, e+e^ reactions, on the one hand, and in central 
or semi-peripheral collisions of A « 200 nuclei, on the other. Comparing in 
detail we formulate a nuclear modification factor for the bulk hadron rapidity 
distributions, 

J.AA ^ dN'^'^/dyjy) inA + A 
y 0.5 Npart dN-h/dy in pp ^ ' 

where Npart < 2^4 is the mean number of "participating nucleons" (which 
undergo at least one inelastic collision with another nucleon) which increases 

with collision centrality. For identical nuclei colliding ^iV^™/^ ~ (-^part ) ^"^^ 

thus 0.5 Npart gives the number of opposing nucleon pairs. R^^ = 1 if each 
such " opposing" pair contributes the same fraction to the total A+ A yield as 
is produced in minimum bias pp at similar ^/s. From Figs. 6 and 8 we infer 
that for I 7y |< 1, R^^ = 1.5 at top RHIC energy, and for the pseudo-rapidity 
integrated total N'^'^ we find R^^ = 1.36, in central Au-|-Au collisions. AA col- 
lisions thus provide for a higher stopping power than pp (which is also reflected 
in the higher rapidity shift {5y) of Fig. 12). The observation that their stop- 
ping power resembles the e+e" inelasticity suggests a substantially reduced 
leading particle effect in central collisions of heavy nuclei. This might not be 
surprising. In a Glauber-view of successive minimum bias nucleon collisions 
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occuring during interpenetration, each participating nucleon is struck v > 3 
times on average, which might saturate the possible inelasticity, removing the 
leading fragment. 

This view naturally leads to the scaling of the total particle production 
in nuclear collisions with N parti as seen clearly in Fig. 8, reminiscent of the 
"wounded nucleon model" [49] but with the scaling factor determined by e+e" 
rather than pp [50] . Overall we conclude from the still rather close similarity 
between nuclear and elementary collisions that the mechanisms of longitudi- 
nal phase space population occur primordially, during interpenetration which 
is over after 0.15 fm/c at RHIC, and after 1.5 fm/c at SPS energy. I.e. it 
is the primordial non-equilibrium pQCD shower evolution that accounts for 
stopping, and its time extent should be a lower limit to the formation time 
To employed in the Bjorken model [34], equation 82. Equilibration at the 
partonic level might begin at t > tq only (the development toward a quark- 
gluon-plasma phase), but the primordial parton redistribution processes set 
the stage for this phase, and control the relaxation time scales involved in 
equilibration [51]. More about this in section 3.5. We infer the existence of a 
saturation scale [52] controlling the total inelasticity: with ever higher reactant 
thickness, proportional to A^^^ , one does not get a total rapidity or energy den- 
sity proportional to A^^^ (the number of " successive binary collisions" ) but to 
^1.08 Qj-^jy jg3j_ jsfote that the lines shown in Fig. 9 (right panel) refer to such a 
saturation theory: the color glass condensate (CGC) model [43] developed by 
McLerran and Venugopulan. The success of these models demonstrates that 
"successive binary baryon scattering" is not an appropriate picture at high 
s/s. One can free the partons from the nucleonic parton density distributions 
only once, and their corresponding transverse areal density sets the stage for 
the ensuing QCD parton shower evolution [52]. Moreover, an additional satu- 
ration effect appears to modify this evolution at high transverse areal parton 
density (see section 3.4). 

3.3 Dependence on system size 

We have discussed above a first attempt toward a variable (Npart) that scales 
the system size dependence in A-l-A collisions. Note that one can vary the 
size either by centrally colliding a sequence of nuclei, Ai + Ai, A2 + A2 etc., 
or by selecting different windows in Npart out of minimum bias collision en- 
sembles obtained for heavy nuclei for which BNL employs ^^"^ Au and CERN 
'^^^Pb. The third alternative, scattering a relatively light projectile, such as 
^"^S, from increasing A nuclear targets, has been employed initially both at the 
AGS and SPS but got disfavored in view of numerous disadvantages, of both 
experimental (the need to measure the entire rapidity distribution, i.e. lab 
momenta from about 0.3-100 GeV/c, with uniform efficiency) and theoretical 
nature (difii'erent density distributions of projectile and target; occurence of 
an"effectiv" center of mass, different for hard and soft collisions, and depend- 
ing on impact parameter). 
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The determination of Npart is of central interest, and thus we need to look 
at technicalities, briefly. The approximate linear scaling with Npart that we 
observed in the total transverse energy and the total charged particle number 
(Figs. 5, 8) is a reflection of the primordial redistribution of partons and 
energy. Whereas all observable properties that refer to the system cvohition 
at later times, which are of interest as potential signals from the equilibrium, 
QCD plasma "matter" phase, have different speciflc dependences on Npart, be 
it suppressions (high px signals, jets, quarkonia production) or enhancements 
(collective hydrodynamic flow, strangeness production). Npart thus emerges 
as a suitable common reference scale. 

Npari captures the number of potentially directly hit nucleons. It is es- 
timated from an eikonal straight trajectory Glauber model as applied to 
the overlap region arising, in dependence of impact parameter b, from the 
superposition along beam direction of the two initial Woods-Saxon density 
distributions of the interacting nuclei. To account for the dilute surfaces of 
these distributions (within which the intersecting nucleons might not find an 
interaction partner) each incident nucleon trajectory gets equipped with a 
transverse radius that represents the total inelastic NN cross section at the 
corresponding ^/s. The formalism is imbedded into a Monte Carlo simulation 
(for detail see [54]) starting from random microscopic nucleon positions within 
the transversely projected initial Woods-Saxon density profiles. Overlapping 
cross sectional tubes of target and projectile nucleons are counted as a partic- 
ipant nucleon pair. Owing to the statistics of nucleon initial position sampling 
each considered impact parameter geometry thus results in a probability dis- 
tribution of derived Npart- Its width a defines the resolution A{b) of impact 
parameter b determination within this scheme via the relation 

^ ,(,)«^%iM (88) 



A{b) ' ' db 

which, at A=200, leads to the expectation to determine b with about 1.5 fm 
resolution [54], by measuring Npart- 

How to measure Npart'^ In fixed target experiments one can calorimetri- 
cally count all particles with beam momentum per nucleon and superimposed 
Fermi momentum distributions of nucleons, i.e. one looks for particles in the 
beam fragmentation domain 2/f,eam±0.5, pt < 0.25GeV/c. These are identified 
as spectator nucleons, and Npart = ^ ~ ^spec- identical nuclear collision 

systems (^N^l°(^ = (Npalf), and thus Npart gets approximated by 2 N^art ■ 
This scheme was employed in the CERN experiments NA49 and WA80, and 
generalized [55] in a way that is illustrated in Fig. 13. 

The top panel shows the minimum bias distribution of total energy registered 
in a forward calorimeter that covers the beam fragment domain in Pb-|-Pb 

collisions at lab. energy of 158 GeV per projectile nucleon, y/s = 17.3 GeV. 
The energy spectrum extends from about 3 TeV which corresponds to about 
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Fig. 13. (a) Energy spectrum of the forward calorimeter in Pb+Pb collisions at 
158^ GeV; (b) impact parameter and fraction of total inelastic cross section related 
to forward energy from the VENUS model [55, 56]. 



20 projectile spectators (indicating a "central" collision), to about 32 TeV 
which is close to the total beam energy and thus corresponds to extremely 
peripheral collisions. Note that the shape of this forward energy spectrum 
is the mirror image of the minimum bias transverse energy distribution of 
Fig. 5, both recorded by NA49. Prom both figures we see that the ideal head- 
on, 6^0 collision can not be selected from these (or any other) data, owing 
to the facts that 6 = carries zero geometrical weight, and that the diffuse 
Woods-Saxon nuclear density profiles lead to a fluctuation of participant nu- 
cleon number at given finite b. Thus the Npart fluctuation at finite weight 
impact parameters overshadows the genuinely small contribution of near zero 
impact parameters. Selecting "central" collisions, either by an on-line trigger 
cut on minimal forward energy or maximal total transverse energy or charged 
particle rapidity density, or by corresponding off-line selection, one thus faces 
a compromise between event statistics and selectivity for impact parameters 
near zero. In the example of Fig. 13 these considerations suggest a cut at 
about 8 TeV which selects the 5% most inelastic events, from among the 
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overall minimum bias distribution, then to be labeled as "central" collisions. 
This selection corresponds to a soft cutoff at 6 < 3 fm. 

The selectivity of this, or of other less stringent cuts on collision centrality 
is then established by comparison to a Glauber or cascade model. The bot- 
tom panel of Fig. 13 employs the VENUS hadron/string cascade model [56] 
which starts from a Monte Carlo position sampling of the nucleons imbed- 
ded in Woods-Saxon nuclear density profiles but (unlike in a Glauber scheme 
with straight trajectory overlap projection) following the cascade of inelastic 
hadron/string multiplication, again by Monte Carlo sampling. It reproduces 
the forward energy data reasonably well and one can thus read off the aver- 
age impact parameter and participant nuclcon number corresponding to any 
desired cut on the percent fraction of the total minimum bias cross section. 
Moreover, it is clear that this procedure can also be based on the total min- 
imum bias transverse energy distribution, Fig. 5, which is the mirror image 
of the forward energy distribution in Fig. 13, or on the total, and even the 
mid-rapidity charged particle density (Fig. 8). The latter method is employed 
by the RHIC experiments STAR and PHENIX. 

How well this machinery works is illustrated in Fig. 14 by RHIC-PHOBOS 

results at ^/s = 200 GeV [42]. The charged particle pseudo-rapidity density 
distributions are shown for central (3-6% highest Nch cut) Cu-|-Cu collisions, 
with (Npart) = 100, and semi-peripheral Au-|-Au collisions selecting the cut 
window (35-40%) such that the same (Npart) emerges. The distributions are 
nearly identical. In extrapolation to Npart = 2 one would expect to find agree- 
ment between min. bias p+p, and "super-peripheral" A-l-A collisions, at least 
at high energy where the nuclear Fermi momentum plays no large role. Fig. 15 
shows that this expectation is correct [57] . As it is technically difficult to select 
^part = 2 from A=200 nuclei colliding, NA49 fragmented the incident SPS 
Pb beam to study ^^C+^^C and '^^Si+'^^ Si collisions [55]. These systems arc 
isospin symmetric, and Fig. 15 thus plots 0.5((7r+) -|- {it~))/ (Nw) including 
p + p where A^^^ = 2 by definition. We see that the pion multiplicity of A-|-A 
collisions interpolates to the p-|-p data point. 

Note that NA49 employs the term "wounded nucleon" number (Nw) to count 

the nucleons that imderwcnt at least one inelastic nucleon-nucleon collision. 
This is what the RHIC experiments (that follow a Glauber model) call Npart 
whereas NA49 reserves this term for nucleons that underwent any inelastic col- 
lision. Thus Nw in Fig. 15 has the same definition as Npart in Figs. 6, 8, 10, 14. 
We see that a smooth increase joins ihep-\-p data, via the light A-|-A central 
collisions, to a saturation setting in with semi-peripheral Pb-|-Pb collisions, 
the overall, relative increase amounting to about 40% (as we saw in Fig. 6). 

There is nothing like an A^^^ increase (the thickness of the reactants) 
observed here, pointing to the saturation mechanism(s) mentioned in the pre- 
vious section, which are seen from Fig. 15 to dampen the initial, fast increase 
once the primordial interaction volume contains about 80 nucleons. In the 
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Fig. 14. Charged hadron pseudo- 
rapidity distributions in Cu+Cu and 
Au-I-Au collisions at a/s = 200 GeV, 
with similar Npart ~ 100 [42]. 



Fig. 15. Charged pion multiplicity 
normalized by Nw vs. centrality in 
p+p, C+C, Si-HSi and Pb-hPb colli- 
sions at = 17.3 GeV [55, 57]. 



Glauber model view of successive collisions (to which we attach only symbol- 
ical significance at high ^/s) this volume corresponds to (i/) « 3, and within 
the terminology of such models we might thus argue, intuitively, that the ini- 
tial geometrical cross section, attached to the nucleon structure function as a 
whole, has disappeared at (v) « 3, all constituent partons being freed. 

3.4 Gluon Saturation in A+A Collisions 

We will now take a closer look at the saturation phenomena of high energy 
QCD scattering, and apply results obtained for deep inelastic electron-proton 
reactions to nuclear collisions, a procedure that relies on a universality of high 
energy scattering. This arises at high y^, and at relatively low momentum 
transfer squared (the condition governing bulk charged particle production 
near mid-rapidity at RHIC, where Feynman x w 0.01 and < 5 GeV"^). 
Universality comes about as the transverse resolution becomes higher and 
higher, with Q^, so that within the small area tested by the collision there 
is no difference whether the partons sampled there belong to the transverse 
gluon and quark density projection of any hadron species, or even of a nucleus. 
And saturation arises once the areal transverse parton density exceeds the 
resolution, leading to interfering QCD sub-amplitudes that do not reflect in 
the total cross section in a manner similar to the mere summation of separate'', 
resolved color charges [51, 52, 53, 43, 58, 59, 60]. 

The ideas of saturation and universality are motivated by HERA deep 
inelastic scattering (DIS) data [61] on the gluon distribution function shown 

in Fig. 16 (left side). The gluon rapidity density, xG{x,Q'^) = ^^^^^ — rises 

* Note that QCD considers interactions only of single charges or charge-anticharge 
pairs. 
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rapidly as a function of decreasing fractional momentum, x, or increasing 
resolution, Q^. The origin of this rise in the gluon density is, ultimately, the 
non-abelian nature of QCD. Due to the intrinsic non-linearity of QCD [59, 60], 
gluon showers generate more gluon showers, producing an avalanche toward 
small X. As a consequence of this exponential growth the spatial density of 
gluons (per unit transverse area per unit rapidity) of any hadron or nucleus 
must increase as x decreases [58]. This follows because the transverse size, 
as seen via the total cross section, rises more slowly toward higher energy 
than the number of gluons. This is illustrated in Fig. 16 (right side). In a 
head-on view of a hadronic projectile more and more partons (mostly gluons) 
appear as x decreases. This picture reflects a representation of the hadron in 
the "infinite momentum frame" where it has a large light-cone longitudinal 
momentum P"*" ^ M. In this frame one can describe the hadron wave function 
as a collection of constituents carrying a fraction p+ = a;P+, < 2; < 1, of the 
total longitudinal momentum [62] ("light cone quantization" method [63]). In 
DIS at large ^/s and one measures the quark distributions dNq/dx at small 
x, deriving from this the gluon distributions xG{x,Q'^) of Fig. 16. 

It is useful [64] to consider the rapidity distribution implied by the parton 
distributions, in this picture. Defining y — yhadron — ln(\/x) as the rapidity 
of the potentially struck parton, the invariant rapidity distribution results as 

dN/dy = X dN/dx = xG(x, Q^). (89) 

At high the measured quark and gluon structure functions are thus simply 
related to the number of partons per unit rapidity, resolved in the hadronic 
wave function. 




Fig. 16. (left) The HERA data for the gluon distribution function as a function of 
fractional momentum x and square momentum transfer [61]. (right) Saturation 
of gluons in a hadron; a head on view as x decreases [64]. 



The above textbook level [63, 64] recapitulation leads, however, to an impor- 
tant application: the dN / dy distribution of constituent partons of a hadron (or 



Global Properties of Nucleus-Nucleus Collisions 



43 



nucleus), determined by the DIS experiments, is similar to the rapidity distri- 
bution of produced particles in hadron-hadron or A+A colHsions as wc expect 
the initial gluon rapidity density to be represented in the finally observed, pro- 
duced hadrons, at high ^/s. Due to the longitudinal boost invariance of the 
rapidity distribution, we can apply the above conclusions to hadron-hadron or 
A-l-A collisions at high ^/s, by replacing the infinite momentum frame hadron 
rapidity by the center of mass frame projectile rapidity, Uproj , while retaining 
the result that the rapidity density of potentially interacting partons grows 
with increasing distance from Uproj like 

= Vproj - y = ln{l/x). (90) 

At RHIC energy, = 200 GeV ., Ay at mid-rapidity thus corresponds to 
X < 10~^ (well into the domain of growing structure function gluon density, 
Fig. 16), and the two intersecting partonic transverse density distributions 
thus attempt to resolve each other given the densely packed situation that 
is depicted in the lower circle of Fig. 16 (right panel). At given (which 
is modest, < 5 GeV"^, for bulk hadron production at mid-rapidity) the 
packing density at mid-rapidity will increase toward higher ^/s as 

^ymidrap ^ ln(^/M), i.e. 1/x « y/^/M (91) 

thus sampling smaller x domains in Fig. 16 according to equation 90. It 
will further increase in proceeding from hadronic to nuclear reaction part- 
ners A-l-A. Will it be in proportion to A^^^l We know from the previous 
sections (3.2 and 3.3) that this is not the case, the data indicating an increase 
with A^-^^. This observation is, in fact caused by the parton saturation effect, 
to which we turn now. 

For given transverse resolution Q"^ and increasing 1/x the parton density 
of Fig. 16 becomes so large that one can not neglect their mutual interactions 
any longer. One expects such interactions to produce " shadowing" , a decrease 
of the scattering cross section relative to incoherent independent scattering 
[59, 60]. As an effect of such shadowed interactions there occurs [64] a satura- 
tion [51, 52, 53, 43, 58, 59, 60, 64] of the cross section at each given Q^, slowing 
the increase with 1/x to become logarithmic once 1/x exceeds a certain crit- 
ical value Xs{Q^). Conversely, for fixed x, saturation occurs for transverse 
momenta below some critical Q^(x), 

Q^W-^^^c^f (92) 

where dN/dy is the a;-dependent gluon density (at y = yproj — ln{l/x)). 

is called the saturation scale. In equation 92 ttR^ is the hadron area (in 
transverse projection), and a^A^c is the color charge squared of a single gluon. 
More intuitively, (x) defines an inversely proportional resolution area Fg (x) 
and at each x we have to choose Fs{x) such that the ratio of total area ttR"^ 
to Fs{x) (the number of resolved areal pixels) equals the number of single 
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gluon charge sources featured by the total hadron area. As a consequence the 
saturation scale Q'^{x) defines a critical arcal resolution, with two different 
types of QCD scattering theory defined, at each x, for > Ql and < Ql, 
respectively [52, 58, 64]. 

As one expects a soft transition between such theories, to occur along the 
transition line implied by Ql{x), the two types of QCD scattering are best 
studied with processes featuring typical well above, or below Ql{x). Jet 
production at ^/s > 200 GeV in pp or AA collisions with typical above 
about 10'^ GeV^, clearly falls into the former class, to be described e.g. by 
perturbative QCD DGLAP evolution of partonic showers [65]. The accronym 
DGLAP refers to the inventors of the perturbative QCD evolution of parton 
scattering with the "runing" strong coupling constant as{Q^), Dokshitzer, 
Gribov, Levine, Altarelli and Parisi. On the other hand, mid-rapidity bulk 
hadron production at the upcoming CERN LHC facility {y/s — lATeV ior pp, 
and 5.5 TeV for A+A), with typical < 5 GeV"^ at x < 10'^, will present 
a clear case for QCD saturation physics, as formulated e.g. in the "Color 
Glass Condensate (CGC)" formalism developed by McLerran, Venugopalan 
and collaborators [43, 58, 64, 66]. This model develops a classical gluon field 
theory for the limiting case of a high areal occupation number density, i.e for 
the conceivable limit of the situation depicted in Fig. 16 (right hand panel) 
where the amalgamating small x gluons would overlap completely, within 
any finite resolution area at modest Q^. Classical field theory captures, by 
construction, the effects of color charge coherence, absent in DGLAP parton 
cascade evolution theories [64]. This model appears to work well already at 
\/s as "low" as at RHIC, as far as small Q"^ bulk charged particle production 
is concerned. We have illustrated this by the CGC model fits [43] to the 
PHOBOS charged particle rapidity distributions, shown in Fig. 9. 

Conversely, QCD processes falling in the transition region between such 
limiting conditions, such that typical w Q'^(x), should present observables 
that are functions of the ratio between the transferred momentum and 
the appropriate saturation scale, expressed by Ql{x). As defines the eff'ec- 
tive transverse sampling area, and Q'^{x) the characteristic areal size at which 
saturation is expected to set in, a characteristic behaviour of cross sections, 
namely that they are universal functions of Q^/Q^, is called "geometric scal- 
ing". The HERA ep scattering data obey this scaling law closely [67], and the 
idea arises to apply the universality principle that we mentioned above: at 
small enough x, all hadrons or nuclei are similar, their specific properties only 
coming in via the appropriate saturation scales Q'^{x, h) or (5^(.t, A). Knowing 
the latter for RHIC conditions we will understand the systematics of charged 
particle production illustrated in the previous section, and thus also be able 
to extrapolate toward LHC conditions in pp and AA collisions. 

All data for the virtual photo-absorption cross section a'^P{x,Q'^) in deep 
inelastic ep scattering with x < 0.01 (which is also the RHIC mid-rapidity 
a;-domain) have been found [67] to lie on a single curve when plotted against 
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Fig. 17. (top) Geometric scaling of the virtual photo-absorption cross section a'"' 
on protons; (middle) cross sections for nuclei normalized according to equation 94; 
(bottom) the ratio of a"'^ to a fit of a~'^ (see [53] for data reference). 



QVQl with 

Qlix) ~ (^)^ IGeV^ (93) 

with A ~ 0.3 and xq — 10"^. This scaling [68] with t ~ Q^/Ql is shown in 
Fig. 17 (top panel) to interpolate all data. A chain of arguments, proposed by 
Armesto, Salgado and Wiedemann [53] connects a fit to these data with photo- 
absorption data for (virtual) photon-A interactions [69] via the geometrical 
scaling ansatz 

TrR\ -kRI ^ ' 

assuming that the scale in the nucleus grows with the ratio of the transverse 
parton densities, raised to the power 1/5 (a free parameter). 

Fig. 17 (middle and bottom panels) shows their fit to the nuclear photo- 
absorbtion data which fixes 6 = 0.79 and ttR^ — 1.57 fm^ (see ref. [53] for 
detail). The essential step in transforming these findings to the case of A+A 
collisions is then taken by the empirical ansatz 
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dy 



(96) 



by which the mid-rapidity parton (ghion) density dN/dy in equation 92 gets 
related to the charged particle mid-rapidity density at y ~ [59, 70], measured 
in nucleus-nucleus collisions. Replacing, further, the total nucleon number 2A 
in a collision of identical nuclei of mass A by the number Npart of participating 
nucleons, the final result is [53] 



where the exponent a = (1 — (5)/35 = 0.089, and Nq = 0.47. The exponent a 

is far smaller than 1/3, a value that represents the thickness of the rcactants, 
and would be our naive guess in a picture of "successive" independent nu- 
cleon participant collisions, whose average number (v) oc {Npartf^Y^^- The 
observational fact (sec Fig. 15) that a < 1/3 for mid-rapidity low bulk 
hadron production in A-|-A collisions illustrates the importance of the QCD 
saturation effect. This is shown [53] in Fig. 18 where equation 97 is applied to 
the RHIC PHOBOS data for mid-rapidity charged particle rapidity density 
per participant pair, in Au-|-Au collisions at y/s = 19.6, 130 and 200 GeV 
[71], also including a prediction for LHC energy. Note that the factorization 
of energy and centrality dependence, implied by the RHIC data [42], is well 
captured by equation 92 and the resulting fits in Fig. 18. Furthermore, the 
steeper slope, predicted for Npart < 60 (not covered by the employed data 
set), interpolates to the corresponding and pp data, at Npart = 2. It resem- 
bles the pattern observed in the NA49 data (Fig. 15) for small Npart collisions 
of light A+A systems, at =17—20 GeV, and may be seen, to reflect the 
onset of QCD saturation. Finally we note that the conclusions of the above, 
partially heuristic approach [53], represented by equations 94 to 97, have been 
backed up by the CGC theory of McLerran and Venugopulan [43, 58, 64] , pre- 
dictions of which we have illustrated in Fig. 9. 

Bulk hadron production in AA collisions at high ^/s can be related, via the 

assumption of imivcrsality of high energy QCD scattering, to the phenomenon 
of geometric scaling first observed in HERA deep inelastic ep cross sections. 
The underlying feature is a QCD saturation effect arising from the diverging 
areal parton density, as confronted with the limited areal resolution Q^, in- 
herent in the considered scattering process. The "saturation scale" Ql{x,A) 
captures the condition that a single partonic charge source within the trans- 
verse partonic density profile can just be resolved by a sufficiently high Q^. 
Bulk hadron production in A-|-A collisions falls below this scale. 

3.5 Transverse phase space: equilibrium and the QGP state 

At RHIC energy, y/s — 200 GeV^, the Au+Au collision rcactants arc longitudi- 
nally contracted discs. At a nuclear radius R « A^^^ fm and Lorentz 7 Ki 100 



1 dN^ 
Npart dy 



{aty^O) = No{^sfNl 



part 



(97) 
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Fig. 18. Saturation model fit [53] applied to RHIC charged hadron multiplicity 
data at mid-rapidity normalized by number of participant pairs, at various energies 
[71]. Also shown is an extrapolation to pp data and a prediction for minimum bias 
Pb-I-Pb collisions at LHC energy, y/s = 5500 GeV . 



their primordial interpenetration phase ends at time tq < 0.15/m/c. This time 
scale is absent in e^e~ annihilation at similar ^/s where tq « 0.1 fm/c marks 
the end of the primordial pQCD partonic shower evolution [72] during which 
the initially created qq pair, of "virtually" Q = ^fsjl each, multiplies in the 
course of the QCD DGLAP evolution in perturbative vacuum, giving rise to 
daughter partons of far lower virtuality, of a few GeV . In A+A collisions this 
shower era should last longer, due to the interpenetrational spread of primor- 
dial collision time. It should be over by about 0.25 fm/c. The shower partons 
in e+e~ annihilation stay localized within back to back cone geometry reflect- 
ing the directions of the primordial quark pair. The eventually observed "jet" 
signal, created by an initial of order IQ'^GeV^, is established by then. Upon 
a slow-down of the dynamical evolution time scale to r « 1 fm/c the shower 
partons fragment further, acquiring transverse momentum and yet lower vir- 
tuality, then to enter a non perturbative QCD phase of color neutralization 
during which hadron-like singlet parton clusters are formed. Their net initial 
pQCD virtuality, in pQCD vacuum, is recast in terms of non-perturbative 
vacuum hadron mass. The evolution ends with on-shell, observed jet-hadrons 
after about 3 fm/c of overall reaction time. 



48 Michael Kliemant, Raghunath Sahoo, Tim Schuster, and Reinhard Stock 

Remarkably, even in this, somehow most elementary process of QCD evo- 
lution, an aspect of equilibrium formation is observed, not in the narrowly 
focussed final dijet momentum topology but in the relative production rates 
of the various created hadronic species. This so-called " hadrochemical" equi- 
librium among the hadronic species is documented in Fig. 19. The hadron 
multiplicities per e+e^ annihilation event at ^/s = 91.2 GeV [4] are con- 
fronted with a Hagedorn [4] canonical statistical Gibbs ensemble prediction 
[73] which reveals that the apparent species equilibrium was fixed at a tem- 
perature of T = 165 MeV , which turns out to be the universal hadronization 
temperature of all elementary and nuclear collisions at high ^/s (Hagedorns 
limiting temperature of the hadronic phase of matter). We shall return to this 
topic in section 4 but note, for now, that reactions with as few as 20 charged 
particles exhibit such statistical equilibrium properties. 
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Fig. 19. Hadron multiplicities in LEP e+e annihilation at ^/s = 91.2 GeV con- 
fronted with the predictions of the canonical statistical hadronization model [73]. 



What happens with parton (and hadron) dynamics in A-f A collisions after 
To? There will not be a QCD evolution in vacuum as the transverse radius 
of the interacting system is large. It may grow to about twice the nuclear 
radius, i.e. to about 15 /m before interactions cease; i.e. the system needs 
about 15 fm/c to decouple. This simple fact is the key to our expectation 
that the expansive evolution of the initial high energy density deposited in a 
cylinder of considerable diameter (about 10 /m-), may create certain equilib- 
rium properties that allow us to treat the contained particles and energy in 
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terms of thermodynamic phases of matter, such as a partonic QGP liquid, or a 
hadronic liquid or gas, etc.. Such that the expansion dynamics makes contact 
to the phase diagram illustrated in Fig. 1. This expectation turns out to be 
justified as we shall describe in section 4. What results for the evolution after 
To in a central A+A collision is sketched in Fig. 2 by means of a schematic 
2-dimensional light cone diagram, which is entered by the two reactant nuclei 
along ±z = t trajectories where z is the beam direction and Lorentz con- 
traction has been taken to an extreme, such that there occurs an idealized 
t = z = interaction "point". Toward positive t the light cone proper time 
profiles of progressing parton-hadron matter evolution are illustrated. The 
first profile illustrated here corresponds to the end of shower formation time 
tq. From our above discussion of the e+e" annihilation process one obtains a 
first estimate, tq > 0.25 fm/c (including interpenetration time of 0.15 fm/c 
at RHIC) which refers to processes of very high > 10'^ GeV^, far above 
the saturation scale Ql discussed in the previous section. The latter scale has 
to be taken into account for low pr hadron production. 

It is the specific resolution scale of a QCD sub-process, as enveloped 
in the overall collision dynamics of two slabs of given transverse partonic 
structure function density, that determines which fraction of the constituent 
partons enters interaction. In the simple case of extremely high processes 
the answer is that all constituents are resolved. However, at modest (domi- 
nating bulk hadron production) the characteristic QCD saturation scale Qs(x) 
gains prominence, defined such that processes with < do not exploit 
the initial transverse parton densities at the level of independent single con- 
stituent color field sources (see equation 92). For such processes the proper 
formation time scale, tq, is of order of the inverse saturation momentum [51], 
1/Qs ~ 0.2 fm/c at ^/s = 200 GeV. The first profile of the time evolution, 
sketched in Fig. 2, should correspond to proper time t = tq = 0.25 fm/c at 
RHIC energy. At top SPS energy, ^/s = 17.3 GeV, we can not refer to such 
detailed QCD considerations. A pragmatic approach suggests to take the in- 
terpenetration time, at 7 « 8.5, for guidance concerning the formation time, 
which thus results as tq ~ 1.5 fm/c. 

In summary of the above considerations we assume that the initial partonic 
color sources, as contained in the structure functions (Fig. 16), are spread out 
in longitudinal phase space after light cone proper time f = tq ~ 0.25 fm/c, at 
top RHIC energy, and after tq ~ 1.4 fm/c at top SPS energy. No significant 
transverse expansion has occured at this early stage, in a central collision 
oi A ^ 200 nuclei with transverse diameter of about 12 fm. The Bjorken 
estimate [34] of initial energy density e (equation 82) refers to exactly this 
condition, after formation time tq. In order to account for the finite longitudi- 
nal source size and interpenetration time, at RHIC, we finally put the average 
To ~ 0.3 fm, at ^/s = 200GeV, indicating the "initialization time" after which 
all partons that have been resolved from the structure functions are engaged 
in shower multiplication. As is apparent from Fig. 2, this time scale is Lorentz 
dilated for partons with a large longitudinal momentum, or rapidity. This 
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means that the slow particles are produced first toward the center of the col- 
lision region, and the fast (large rapidity) particles arc produced later, away 
from the collision region. This Bjorken "inside-out" correlation [34] between 
coordinate- and momentum-space is similar to the Hubble expansion pattern 
in cosmology: more distant galaxies have higher outward velocities. Analo- 
gously, the matter created in A-|-A collisions at high ^/s is born expanding, 
however with the difference that the Hubble flow is initially one dimensional 
along the collision axis. This pattern will continue, at ^/s = 200 GeV , until 
the system begins to feel the effects of finite size in the transverse direction 
which will occur at some time to in the vicinity of 1 fm/c. However, the tight 
correlation between position and momentum initially imprinted on the system 
will survive all further expansive evolution of the initial " firetube" , and is well 
recovered in the expansion pattern of the finally released hadrons of modest 
Pt as wc shall show when discussing radial flow (sec section 3.6). 

In order to proceed to a more quantitative description of the primordial 
dynamics (that occurs onward from tq for as long the time period of pre- 
dominantly longitudinal expansion might extend) wc return to the Bjorken 
estimate of energy density, corresponding to this picture [34], as implied by 
equation 82, which we now recast as 

e={^^){El){.R\to)-' (98) 

where the first term is the (average) total hadron multiplicity per unit rapid- 
ity which , multiplied with the average hadron transverse energy, equals the 
total transverse energy recorded in the calorimetric study shown in Fig. 5, as 
employed in equation 82. The quantity Ra is, strictly speaking, not the radius 
parameter of the spherical Woods-Saxon nuclear density profile but the rms 
of the reactant overlap profiles as projected onto the transverse plane (and 
thus slightly smaller than Ra ~ A'^/^ fm). Employing A^^^ here (as is done 
throughout) leads to a conservative estimate of e, a minor concern. However, 
the basic assumption in equation 98 is to identify the primordial transverse en- 
ergy "radiation", of an interactional cylindric source of radius Ra and length 
^0 (where Tq < to < 1 fm/c, not Lorentz dilated at midrapidity) , with the 
finally emerging bulk hadronic transverse energy. We justify this assumption 
by the two observations, made above, that 

1. the bulk hadron multiplicity density per unit rapidity resembles the 
parton density, primordially released at saturation scale tq (Figs. 9, 18) 
at y/s = 200 GeV, and that 

2. the global emission pattern of bulk hadrons (in rapidity and pt) closely 
reflects the initial correlation between coordinate and momentum space, 
characteristic of a primordial period of a predominantly longitudinal ex- 
pansion, as implied in the Bjorken model. 

Both these observations are surprising, at first sight. The Bjorken model was 
conceived for elementary hadron collisions where the expansion proceeds into 
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vacuum, i.e. directly toward observation. Fig. 2 proposes that, to the contrary, 
primordially produced partons have to transform through further, successive 
stages of partonic and hadronic matter, at decreasing but still substantial en- 
ergy density, in central A-|-A collisions. The very fact of high energy density, 
with implied short mean free path of the constituent particles, invites a hydro- 
dynamic description of the expansive evolution. With initial conditions fixed 
between tq and to, an ensuing 3-dimensional hydrodynamic expansion would 
preserve the primordial Bjorkcn-typc correlation between position and mo- 
mentum space, up to lower density conditions and, thus, close to emission of 
the eventually observed hadrons. We thus feel justified to employ equation 82 
or 98 for the initial conditions at RHIC, obtaining [51, 73] 

6 GeV/fm^ < e < 20 GeV/fm^ (99) 

for the interval 0.3 /m/c < to < 1 fm/c, in central Au-|-Au collisions at y w 
and ^/s = 200 GeV. The energy density at top SPS energy, ^/s = 17.3 GeV, 
can similarly be estimated [31, 32] to amount to about 3 GeV/fm^ at a to 
of 1 fm/c but we can not identify conditions at tq < to in this case as the 
mere interpenetration of two Pb nuclei takes 1.4 fm/c. Thus the commonly 
accepted to — 1 fm/c may lead to a high estimate. An application of the 
parton-hadron transport model of Ellis and Geiger [74, 75] to this collision 
finds e = 3.3 GeV/ fm^ &t t = 1 fm/c. A primordial energy density of about 
3 GeV/ fm^ is twenty times po » 0.15 GeV/ fm^, the average energy density 
of ground state nuclear matter, and it also exceeds, by far, the critical QCD 
energy density, of 0.6 < Cc < 1 GeV/ fm^ according to lattice QCD [38]. The 
initial dynamics thus clearly proceeds in a deconfined QCD system also at top 
SPS energy, and similarly so with strikingly higher energy density, at RHIC, 
where time scales below 1 fm/c can be resolved. 

However, in order now to clarify the key question as to whether, and when 
conditions of partonic dynamical equilibrium may arise under such initial 
conditions, we need estimates both of the proper relaxation time scale (which 
will, obviously, depend on energy density and related collision frequency), 
and of the expansion time scale as governed by the overall evolution of the 
collision volume. Only if T{relax.) < T{expans.) one may conclude that the 
"deconfined partonic system" can be identified with a "deconfined QGP state 
of QCD matter" as described e.g. by lattice QCD, and implied in the phase 
diagram of QCD matter suggested in Fig. 1. 

For guidance concerning the overall time-order of the system evolution 
wc consider information [76] obtained from Bosc-Einstein correlation anal- 
ysis of pion pair emission in momentum space. Note that pions should be 
emitted at any stage of the evolution, after formation time, from the surface 
regions of the evolving "fire-tube". Bulk emission of pions occurs, of course, 
after hadronization (the latest stages illustrated in the evolution sketch given 
in Fig. 2). The dynamical pion source expansion models by Heinz [77] and 
Sinyukov [78] elaborate a Gaussian emission time profile, with mean r/ (the 
decoupling time) and width At (the duration of emission). 
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Fig. 20. Time profile of pion decoupling rate from the fireball in a central Pb+Pb 
collision, with r = the end of the formation phase. Bose-Einstein correlation of 
7r~7r~pairs yields an average Gaussian decoupling profile with r/ = 8 fm/c and 
duration of emission parameter At = 4 fm/c [76, 77]. 



Fig. 20 shows an application of this analysis to central Pb+Pb coUision 
negative pion pair correlation data obtained by NA49 at top SPS energy, 
•^s = 17.3GeV [79], whore rj « 8 fm/c and At « 4 fm/c (note that r = in 
Fig. 20 corresponds, not to interaction time t = but to t ~ lAfm/c, the end 
of the interpenetration phase) . We see, first of all, that the overall dynamical 
evolution of a central Pb+Pb collision at = 17.3 GeV is ending at about 
15 fm/c; the proper time defines the position of the last, decoupling profile 
illustrated in Fig. 2, for the SPS collisions considered here. While the details 
of Fig. 20 will turn out to be relevant to our later discussion of hadronization 
(section 4) and hadronic expansion, we are concerned here with the aver- 
age proper time at which the partonic phase ends. After consideration of the 
duration widths of these latter expansion phases [75, 76] one arrives at an es- 
timate for the average time, spent before hadronization, of Z\t = 3 — 4 fm/c, 
again in agreement with the parton cascade model mentioned above [75]. 
This model also leads to the conclusion that parton thermal equilibrium is, 
at least, closely approached locally in these central Pb+Pb collisions as far as 
mid-rapidity hadron production is concerned (at forward-backward rapidity 
the cascade re-scattering processes do not suffice, however). 

This finding agrees with earlier predictions of Treiax = 1 — 2 fm/c at top 
SPS energy [80]. However we note that all such calculations employ pertur- 
bative QCD methods, implying the paradoxical consequence that equilibrium 
is closely approached toward the end of the partonic phase, at such low ^/s, 
i.e. in a QGP state at about T = 200 MeV which is, by definition, of non- 
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perturbative nature. We shall return to the question of partonic equilibrium 
attainment at SPS energy in the discussion of the hadronization process in 
nuclear collisions (section 4). 

Equilibrium conditions should set in earlier at top RHIC energy. As 
transverse partonic expansion should set in after the proper time interval 
0.3 /m/c <to<l fm/c (which is now resolved by the early dynamics, unlike 
at top SPS energy), we take guidance from the Bjorken estimate of primordial 
energy density which is based on transverse energy production data. Conser- 
vatively interpreting the result in equation 99 we conclude that e is about four 
times higher than at ^/s = 17.3 GeF in the above proper time interval. As the 
binary partonic collision frequency scales with the square of with the square 
of the density p (related to the energy density e via the relation e = {E) p = 
Tp), and is inversely proportional to the relaxation time Trelax we expect 

Trela. CX {1/ pf « {T / ef (100) 

which implies that Treiax{RHIC) ~ 0.25 Treiax{SPS) ~ 0.5 fm/c if we em- 
ploy the estimate T{RHIC) = 2T{SPS). This crude estimate is, however, 
confirmed by the parton transport model of Molar and Gyulassy [81]. 

Partonic equilibration at i/i = 200 GeV should thus set in at a time scale 
commensurate to the (slightly smaller) formation time; scale, at which the to 
be participant partons are resolved from the initial nucleon structure functions 
and enter shower multiplication. Extrapolating to the conditions expected at 
LHC energy (y^ = 5.5 TeV for A+A collisions), where the initial parton 
density of the structure functions in Fig. 16 is even higher {x « 10~^ at mid- 
rapidity), and so is the initial energy density, we may expect conditions at 
which the resolved partons arc almost " born into equilibrium" . 

Early dynamical local equilibrium at RHIC is required to understand 
the observations concerning elliptic flow. This term refers to a collective 
anisotropic azimuthal emission pattern of bulk hadrons in semi-peripheral 
collisions, a hydrodynamical phenomenon that originates from the initial ge- 
ometrical non-isotropy of the primordial interaction zone [82, 83]. A detailed 
hydrodynamic model analysis of the corresponding elliptic flow signal at RHIC 
[84] leads to the conclusion that local equilibrium (a prerequisite to the hy- 
drodynamic description) sets in at to « 0.5 /m/c. This conclusion agrees with 
the estimate via equation 100 above, based on Bjorken energy density and 
corresponding parton collisions frequency. 

We note that the concept of a hydrodynamic evolution appears to be, al- 
most necessarily ingrained in the physics of a system born into (Hubblc-type) 
expansion, with a primordial correlation between coordinate and momentum 
space, and at extreme initial parton density at which the partonic mean free 
path length A is close to the overall spatial resolution resulting from the sat- 
uration scale, i.e. A « l/Qs- 

The above considerations suggest that a quark-gluon plasma state should 
be created early in the c;xpansion (lynaniic;s at \/s = 200 GeV, at about 
T = 300 Mey, that expands hydrodynamically until hadronization is reached, 
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at T w 165 — 170 MeV. Its manifestations will be considered in section 4. At 
the lower SPS energy, up to 17.3 GeV, we can conclude, with some caution, 
that a deconfined hadronie matter system should exist at T fa 200 MeV, in 
the closer vicinity of the hadronization transition. It may closely resemble the 
QGP state of lattice QCD, near Tc- 



3.6 Bulk hadron transverse spectra and radial expansion flow 

In this section wc analyze hulk hadron transverse momentum spectra ob- 
tained at SPS and RHIC energy, confronting the data with predictions of the 
hydrodynamical model of collective expansion matter flow that we have sug- 
gested in the previous section, to arise, almost necessarily, from tlic^ primordial 
Hubble-typc coupling between coordinate and momentum space that prevails 
at the onset of the dynamical evolution in A+A collisions at high ^/s. As 
all hadronie transverse momentum spectra initially follow an approximately 
exponential fall-off (see below) the bulk hadronie output is represented by 
thermal transverse spectra at pT < 2 GeV/c. 

Furthermore we shall focus here on mid-rapidity production in near central 
A-|-A collisions, because hydrodynamic models refer to an initialization period 
characterized by Bjorken-type longitudinal boost invariance, which we have 
seen in Figs. 9 and 11 to be restricted to a relatively narrow interval centered at 
mid-rapidity. Central collisions are selected to exploit the azimuthal symmetry 
of emission, in an ideal impact parameter b geometry. We thus select 
the predominant, relevant hydrodynamic "radial flow" expansion mode, from 
among other, azimuthaly oriented (directed) flow patterns that arise once this 
cylindrical symmetry (with respect to the beam direction) is broken in finite 
impact parameter geometries. 

In order to define, quantitatively, the flow phenomena mentioned above, 
we rewrite the invariant cross section for production of hadron species i in 
terms of transverse momentum, rapidity, impact parameter b and azimuthal 
emission angle ipp (relative to the reaction plane), 

dN,{b) 1 dNM r. „ w ,x „ w ,x ^ n 

PTapTdydipp 2it pTdprdy ■' 

(101) 

where we have expanded the dependence on into a Fourier series. Due to 
reflection symmetry with respect to the reaction plane in collisions of identical 
nuclei, only cosine terms appear. Restricting to mid-rapidity production all 
odd harmonics vanish, in particular the " directed flow" coefficient v\ , and we 
have dropped the y-dependence in the flow coefficients v\ and V2- The latter 
quantifies the amount of "elliptic flow" as discussed above. In the following, we 
will restrict to central collisions which we shall idealize as near-zero impact 
parameter processes governed by cylinder symmetry, whence all azimuthal 
dependence (expressed by the v\, v\,... terms) vanishes, and the invariant 
cross section reduces to the first term in equation 101, which by definition 
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also corresponds to all measurements in which the orientation of the reaction 
plane is not observed. 

Typical transverse momentum spectra of the latter type are shown in 
Fig. 21, for charged hadron production in Au+Au collisions at \fs = 200GeV, 
exhibiting mid-rapidity data at various collision centralities [85]. We observe 
a clear-cut transition, from bulk hadron emission at px < 2 GeV/c featur- 
ing a near-exponential cross section (i.e. a thermal spectrum), to a high pr 
power-law spectral pattern. Within the context of our previous discussion 
(section 3.4) we tentatively identify the low px region with the QCD physics 
near saturation scale. Hadron production at px — > 10 GeV/c should, on the 
other hand, be the consequence of primordial leading parton fragmentation 
originating from "hard", high perturbative QCD processes. 

We thus identify bulk hadron production at low px as the emergence of the 
initial parton saturation conditions that give rise to high energy density and 
small equilibration time scale, leading to a hydrodynamical bulk matter ex- 
pansion evolution. Conversely, the initially produced hard partons, from high 

processes, arc not thcrmalizcd into the bulk but traverse it, as tracers, 
while being attenuated by medium-induced rescattering and gluon radiation, 
the combined effects being reflected in the high px inclusive hadron yield, and 
in jet correlations of hadron emission. We can not treat the latter physics ob- 
servables in detail here, but will remain in the field of low px physics, related 
to hydrodynamical expansion modes, focusing on radially symmetric expan- 
sion. 



In order to infer from the spectral shapes of the hadronic species about the 
expansion mechanism, we first transform to the transverse mass variable, 

'm-T = (Px + rn?y^^, via 

I = ! (102) 

2tt pxdpxdy 27r mxdmxdy 

because it has been shown in p-|-p collisions [86] near RHIC energy that the 
m-x distributions of various hadronic species exhibit a universal pattern ("my 
scaling") at low mx'- 

with a universal inverse slope parameter T and a species dependent normal- 
ization factor A. Hagedorn showed [87] that this scaling is characteristic of 
an adiabatic expansion of a fireball at temperature T. We recall that, on the 
other hand, an ideal hydrodynamical expansion is isentropic. 

Fig. 22 shows the ^/s dependence of the average transverse kinetic energy 
(rriy) — m* for pions, kaons and protons observed at mid-rapidity in central 
Au+Au/Pb+Pb collisions [44]. Similarly, the inverse slope parameter T re- 
sulting from a fit of equation 103 to and K~ transverse mass spectra (at 
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Fig. 21. Transverse momentum spectra of charged hadrons in Au+Au collisions 
at ^/s — 200 GeV, in dependence of collision centrality [85] (offset as indicated), 
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Pt < 2 GeV/c) is shown in Fig. 23, both for nuclear and p+p coUisions [88]. 
We see, first of all, that scaling does not apply in A+A collisions, and that 
the kaon inverse slope parameter, T « 230 MeV over the SPS energy regime, 
can not be identified with the fireball temperature at hadron formation which 
is Th ~ 165 MeV from Fig. 1. The latter is seen, however, to be well repre- 
sented by the p+p spectral data exhibited in the left panel of Fig. 23. There 
is, thus, not only thermal energy present in A-l-A transverse expansion, but 
also hydrodynamical radial flow. 

We note that the indications in Figs. 22 and 23, of a plateau in both (ttit) 
and T, extending over the domain of SPS energies, 6 < -^/s < 17 GeV, have 
not yet been explained by any fundamental expansive evolution model, in- 
cluding hydrodynamics. Within the framework of the latter model, this is a 
consequence of the initialization problem [89] which requires a detailed model- 
ing, both of primordial energy density vs. equilibration time scale, and of the 
appropriate partonic matter equation of state (EOS) which relates expansion 
pressure to energy density. At top RHIC energy, this initialization of hydro- 
flow occurs, both, at a time scale s:! 0.5 fm/c which is far smaller than the 
time scale of eventual bulk hadronization (<: « 3 fm/c), and at a primordial 
energy density far in excess of the critical QCD confinement density. After 
initialization, the partonic plasma phase thus dominates the overall expansive 
evolution, over a time interval far exceeding the formation and relaxation time 
scale. 

Thus, at RHIC energy, parton transport [81] and relativistic hydrodynamic 
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Fig. 23. The inverse slope parameter T of equation 103 for K and transverse 
mass spectra at pr < 2 GeV/c and mid-rapidity in central A-f A, and in minimum 
bias p-f-p collisions [88]. 



[84, 89] models establish a well developed expansion mode that survives the 
subsequent stages of hadronization and hadronic expansion. This is reflected in 
their success in describing elliptic flow. On the other hand, the hydrodynami- 
cal model far overestimates elliptic flow at SPS energy [89] at which, as we have 
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shown in sections. 5, the initiahzation period may be not well separated from 
the confinement (hadronization) stage. Thus, whereas the expansion evolution 
at = 200 GeV (occuring at near-zero baryo-chemical potential in Fig. 1) 
"races" across the parton-hadron phase boundary with fully established flow 
patterns, near fis — where lattice QCD predicts the phase transformation 
to be merely a soft cross-over [18], the dynamics at -^i = 10 — 20 GeV may 
originate from only slightly above, or even at the phase boundary, thus sam- 
pling the domain 200 < /is < 500 MeV where the equation of state might 
exhibit a " softest point" [89] . The hydrodynamic model thus faces formidable 
uncertainties regarding initialization at SPS energy. 

The plateau in Figs. 22, 23 may be the consequence of the fact that not 
much flow is generated in, or transmitted from the partonic phase, at SPS 
energies, because it is initialized close to the phase boundary [88] where the 
expected critical point [11, 12] (Fig. 1), and the corresponding adjacent first 
order phase transition might focus [90] or stall [89] the expansion trajectory, 
such that the observed radial flow stems almost exclusively from the hadronic 
expansion phase. The SPS plateau, which we shall subsequently encounter in 
other bulk hadron variables (elliptic flow, HBT radii) might thus emerge as a 
consequence of the critical point or, in general, of the flatness of the parton- 
hadron coexistence line. RHIC dynamics, on the other hand, originates from 
far above this line. 

Hadronic expansion is known to proceed isentropically [91]: commensu- 
rate to expansive volume increase the momentum space volume must decrease, 
from a random isotropic thermal distribution to a restricted momentum orien- 
tation preferentially perpendicular to the flreball surface, i.e. radial. The initial 
thermal energy, implied by the hadron formation temperature Th = l65MeV, 
will thus fall down to a residual Tp at hadronic decoupling from the flow field 
("thermal freeze-out") plus a radial transverse kinetic energy term rui 
where is the mass of the considered hadron species and (Pt) the average 
radial velocity. We thus expect [92] for the slope of equation 103: 



the latter expression valid at pr larger than hadron mass scale (T then is the 
"blue-shifted temperature" at decoupling [93] and (vt) the average transverse 
velocity) . The assumption that radial flow mostly originates from the hadronic 
expansion phase is underlined by the proportionality of flow energy to hadron 
mass (equation 104). 

Fig. 24 illustrates this proportionality, by a recent compilation [92] of RHIC 
results for central Au-I-Au collisions at i/i = 200 GeV, and SPS results for 

central Pb+Pb collisions at top SPS energy, = 17.3 GeV. At the latter 
energy the slope parameter of the meson is seen to be close to that of the 



T = TF + mi {Prf ,Pt<2 GeV/c 



(104) 



and 




(105) 
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similar mass baryons p and A, emphasizing the occurence of rrii scaling as 
opposed to valence quark number scaling that we will encounter in RHIC 
elliptic flow data [83]. As is obvious from Fig. 24 the multi-strange hyperons 
and charmonia exhibit a slope saturation which is usually explained [92] as 
a consequence of their small total cross sections of rescattering from other 
hadrons, leading to an early decoupling from the bulk hadron radial flow 
field, such that {Pt) ^ < {Pt)^- 

According to our observations with equation 104 a hydrodynamical ansatz 
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Fig. 25. Hadron decoupling temperature T/, and average radial flow velocity {/3t) 
extracted from blast wave model (see equation 106) fits of rriT spectra vs. ^/s [44]. 



for the transverse mass spectrum of hadrons should thus contain the variables 
"true temperature" Tp at decoupling from the flow field, and its average 
velocity {Pt), common to all hadrons. This is the case for the blast wave model 
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[93] developed as an approximation to the full hydrodynamic formalism [89] , 
assuming a common decoupling or "frcczc-out" from flow, for all hadronic 
species, and a boost-invariant longitudinal expansion: 



where p = tanh~^j3T- In an extended version of this model a function is 
included that describes the radial profile of the transverse velocity field, 
I^Tir) = (3^°'^r/R, instead of employing a fixed Pt at decoupHng [94]. Fig. 25 
shows [44] the resulting energy dependence of Tp and (/3t), for the same set 
of data as implied already in Figs. 22 and 23. The "true" decoupling temper- 
ature rises steeply at the AGS and less so at SPS energy (as does (/3t)), to a 
value of about 95 MeV at top SPS energy, which is considerably lower than 
the chemical freeze-out temperature, Th = 165 MeV, at which the hadronic 
species relative yield composition of the hadronic phase becomes stationary 
(see section 4, and Fig. 1). Chemical decoupling thus occurs early, near the 
parton-hadron phase boundary, whereas hadronic radial flow ceases after sig- 
nificant further expansion and cooling, whence the surface radial velocity (its 
average value given by (/3t) in Fig. 25) approaches /3t ~ 0.65. Both data sets 
again exhibit an indication of saturation, over the interval toward top SPS 
energy: the SPS plateau. This supports our above conjecture that radial flow 
is, predominantly, a consequence of isentropic bulk hadronic expansion in this 
energy domain, which sets in at Th. At RHIC energy, both parameters exhibit 
a further rise, suggesting that primordial partonic flow begins to contribute 
significantly to radial flow. 

In summary we conclude that hadronic hadro-chemical freeze-out occurs 
near the parton-hadron boundary line of Fig. 1, at high ^/s. This fixes the 
hadronic species equilibrium distribution at T = Th. Onward from hadronic 
freeze-out, the expanding system is out of equilibrium as far as its hadronic 
species composition is concerned. However the momentum space distributions 
are subject to further change during the hadronic expansion phase, predomi- 
nantly due to elastic scattering and development of radial flow. This process 
ends in final decoupling from strong interaction, at Tf < Th- Analysis of j)t 
spectra exhibits a collective radial flow pattern at Tf, as characterized by a 
radial flow velocity /3t- 
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4 Hadronization and hadronic freeze-out in A+A 
collisions 

Within the course of the global expansion of the primordial reaction volume 
the local flow "cells" will hit the parton-hadron phase boundary as their en- 
ergy density approaches Ccrit ~ 1 GeV/ fw?. Hadronization will thus occur, 
not at an instant over the entire interaction volume, but within a finite overall 
time interval [75] that results from the spread of proper time at which indi- 
vidual cells, or coherent clusters of such cells (as developed during expansion) 
arrive at the phase boundary. However, irrespective of such a local-temporal 
occurence, the hadronization process (which is governed by non perturbative 
QCD at the low corresponding to bulk hadronization) universally results 
in a novel, global equilibrium property that concerns the relative abundance of 
produced hadrons and resonances. This so-called " hadrochemical equilibrium 
state" is directly observable, in contrast to the stages of primordial parton 
equilibration that are only indirectly assessed, via dynamical model studies. 

This equilibrium population of species occurs both in elementary and nu- 
clear collisions [95]. We have seen in Fig. 19 a first illustration, by e+e~ anni- 
hilation data at ^/s = 91.2 GeV LEP energy, that are well reproduced by the 
partition functions of the statistical hadronization model (SHM) in its canon- 
ical form [73]. The derived hadronization temperature, Th = 165 MeV, turns 
out to be universal to all elementary and nuclear collision processes at ^/s > 
20 GeV, and it agrees with the limiting temperature predicted by Hagedorn 
[4] to occur in any multi-hadronic equilibrium system once the energy density 
approaches about Q.QGeV/ fm^ . Thus, the upper limit of hadronic equilibrium 
density corresponds, closely, to the lower limit, ecrit — 0.6 — 1.0 GeV/ fm^ of 
partonic equilibrium matter, according to lattice QCD [38]. In elementary col- 
lisions only about 20 partons or hadrons participate: there should be no chance 
to approach thermodynamic equilibrium of species by rescattering cascades, 
neither in the partonic nor in the hadronic phase. The fact that, nevertheless, 
the hadron formation temperature Th coincides with the Hagedorn limiting 
temperature and with the QCD confinement temperature, is a consequence 
of the non-perturbative QCD hadronization process itself [74], which "gives 
birth" to hadrons/resonances in canonical equilibrium, at high ^/s, as we shall 
see below. This process also governs A+A collisions but, as it occurs here un- 
der conditions of high energy density extended over considerable volume, the 
SHM description now requires a grand canonical ensemble, with important 
consequences for production of strange hadrons (strangeness enhancement). 

The grand canonical order of hadron/resonance production in central A+A 
collisions, and its characteristic strangeness enhancement shows that a state 
of extended matter that is quantum mechanically coherent must exist at 
hadronization [76, 77, 95]. Whether or not it also reflects partonic equilib- 
rium properties (including flavor equilibrium), that would allow us to claim 
the direct observation of a quark gliion plasma state near T,., can not be 
decided on the basis of this observation alone, as the hadronization process 
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somehow generates, by itself, the observed hadronic equiUbrium. This conclu- 
sion, however, is still the subject of controversy [95] . 
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Fig. 26. Total hadron multiplicities in central Pb+Pb collisions at ^/s = 17.3 GeV 
[88] versus prediction of the grand canonical statistical hadronization model [21]. 

Two typical examples of grand canonical SHM application are illustrated in 
Figs. 26 and 27, the first showing total hadron multiplicities in central Pb+Pb 
collisions at ^ = 17.3 GeV by NA49 [88] confronted with SHM predictions 
by Becattini et al. [21]. This plot is similar to Fig. 19 in which e+e~ anni- 
hilation to hadrons is confronted with a SHM prediction derived from the 
canonical ensemble [73]. Central Au+Au collision data at ^/s = 200 GeV 
from several RHIC experiments are compared to grand canonical model pre- 
dictions by Braun-Munzingor ct al. [96] in Fig. 27. The key model parameters, 
Th and the baryo-chemical potential result as 159 MeV {160 MeV), and 
247 MeV (20 MeV) at y/s = 17.3 (200) GeV, respectively. The universahty of 
the hadronization temperature is obvious from comparison of the present val- 
ues with the results of the canonical procedure employed in e+e~ annihilation 
to hadrons a.t y/s = 91.2 GeV (Fig. 19), and in canonical SHM fits [97] to p+p 
collision data at ^ = 27 A GeV where Th = 159 and 169 MeV, respectively. 

Figs. 26 and 27 illustrate two different approaches employed in grand canonical 
SHM application, the former addressing the values of the hadronic multiplici- 
ties as obtained in approximate full 47r acceptance (within limitations implied 
by detector performance), the latter employing a set of multiplicity ratios ob- 
tained in the vicinity of mid-rapidity as invited, at RHIC energy, by the limited 
acceptance of the STAR and PHENIX experiments. The latter approach is ap- 
propriate, clearly, in the limit of boost-invariant rapidity distributions where 
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Fig. 27. Hadron multiplicity ratios at mid-rapidity in central Au-f An collisions at 
y/s = 200 GeV from RHIC experiments STAR, PHENIX and BRAHMS, compared 
to predictions of the grand canonical statistical model [96]. 



hadron production ratios would not depend on the choice of the observational 
rapidity interval. We have shown in section 3.2 that such conditions do, in 
fact, set in at top RHIC energy, as referred to in Fig. 27. However, at low ^/s 
the y-distributions are far from boost-invariant, and the total rapidity gap Ay 
may become comparable, in the extreme case, to the natural rapidity widths 
of hadrons emitted in the idealized situation of a single, isotropically decaying 
fireball positioned at mid-rapidity. Its rapidity spectra, equation 86, resemble 
Gaussians with widths Fi w 2.35 [T/rriiY/'^ for hadron masses nii. Clearly, 
the particle ratios (dNi/dy)/ {dNj/dy) then depend strongly on the position 
of the rapidity interval dy: away from y = heavy hadrons will be strongly 
suppressed, and particle yields in narrow rapidity intervals are useless for a 
statistical model analysis unless it is known a priori that the radiator is a 
single stationary spherical fireball [98]. This is not the case toward top SPS 
energy (see Fig. 12), due to significant primordial longitudinal expansion of 
the hadron emitting source. Given such conditions, the total multiplicity per 
collision event (the invariant yield divided by the total overall inelastic cross 
section) should be employed in the SHM analysis, as is exemplified in Fig. 26. 

4.1 Hadronic freeze-out from expansion flow 

The hadronic multiplicities result from integration of the invariant triple dif- 
ferential cross section over px and y. Instrumental, experiment-specific condi- 
tions tend to result in incomplete px and/or y acceptances. It is important to 
ascertain that the effects of hydrodynamic transverse and longitudinal flow do 
not blast a significant part of the total hadron yield to outside the acceptance. 
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and that they, more generally, do not change the relative hadron yield com- 
position, thus basically affecting the SHM analysis. To see that hadronization 
incorporates only the internal energy in the co-moving frame [98], we first 
assume that hadrochemical freeze-out occurs on a sharp hypersurface S, and 
write the total yield of particle species i as 

^i=J^J^ P'^d^'^^ix) Mx,p) = Jj'a^{x)j^{x) (107) 

where (fia is the outward normal vector on the surface, and 

Jii^) =9vj dV0(/)(5(p2 - mf) p^(exp [p • u{x) - n,]/T ± l)-^ (108) 

is the grand canonical number current density of species z, /i^ the chemical 
potential, u{x) the local flow velocity, and gi the degeneracy factor. In thermal 
equilibrium it is given by 

j1{x) = pi{x)u^{x) with 

Pi{x) = u^{x)j'i{x) = j dV0(p°)^(/ - ml) p ■ u{x) Mp ■ u{x); T; pn) 

= j rfV fi{Ep^;T,pi) = Pi{T,iii). (109) 

Here Epi is the energy in the local rest frame at point x. The total particle 
yield of species i is therefore 

Ni = pi{T,pi) jj''a^{x)u^{x)=Pi{T,pi)Vs{u^) (110) 

where only the total comoving volume Vs of the freeze-out hypersurface S 
depends on the flow profile . V is thus a common total volume factor at 
hadronization (to be determined separately), and the flow pattern drops out 
from the yield distribution over species in 47r acceptance [98]. For nuclear 
collisions at SPS energies and below one thus should perform a SHM analysis 
of the total, 47r-integrated hadronic multiplicities, as was done in Fig. 26. 

We note that the derivation above illustrates the termination problem of 
the hydrodynamic description of A-|-A collisions, the validity of which depends 
on conditions of a short mean free path, A < 1 fm. A precise argumentation 
suggests that two different free paths arc relevant here, concerning hadron oc- 
cupation number and hadron spectral freeze-out, respectively. As hadrochemi- 
cal freeze-out occurs in the immediate vicinity of T,; (and Th ~ 160 — 165Mel^ 
from Figs. 26, 27), the hadron species distribution stays constant through- 
out the ensuing hadronic phase, i.e the "chemical" mean free path abruptly 
becomes infinite at Th, whereas elastic and resonant rescattering may well 
extend far into the hadronic phase, and so does collective pressure and flow. 
In fact we have seen in section 3.6 that the decoupling from flow occurs at Tp 
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as low as 90 — 100 MeV (Fig. 25). Thus the hydrodynamic evolution of high 

collisions has to be, somehow artificially, stopped at the partoii-hadron 
boundary in order to get the correct hadron multiplicities Ni, of equations 107 
to 110, which then stay frozen-out during the subsequent hadronic expansion. 

The equations (107-110) demonstrate the application of the Coopcr-Fryc 
prescription [99] for termination of the hydrodynamic evolution. The hyper- 
surface S describes the space-time location at which individual flow cells 
arrive at the freeze-out conditions, e = ec and T ^ Tc, oi hadronization. At 
this point, the resulting hadron/resonance spectra (for species i) are then 
given by the Cooper-Prye formula 

^§^ = Z^^ = 7#l3 / Mp-u{x),x)p.dM^), (111) 

where p'^fid^cr^ is the local flux of particle i with momentum p through the 

surface S. For the phase space distribution / in this formula one takes the 
local equilibrium distribution at hadronic species freeze-out from the grand 
canonical SHM 

h{E,x) = [exp{{E, - ^^^{x))/T) ± 1]"! (112) 

boosted with the local flow velocity u^{x) to the global reference frame by the 
substitution E ^ p ■ u{x). Fixing T = Tc (taken e.g. from lattice QCD) the 
hadron multiplicities Ni then follow from equation 110, and one compares to 
experiment, as in Figs. 26, 27. In order now to follow the further evolution, 
throughout the hadronic rescattering phase, and to finally compare predic- 
tions of equation 111 to the observed flow data as represented by the various 
Fourier-terms of equation 101 one has to re-initialize (with hadronic EOS) the 
expansion from E{Tc = 165 MeV) untfl final decoupling [89], at T « lOOMel^, 
thus describing e.g. radial and elliptic flow. 

Alternatively, one might end the hydrodynamic description at T = Tc and 
match the thus obtained phase space distribution of equation 111 to a mi- 
croscopic hadron transport model of the hadronic expansion phase [84, 100]. 
This procedure is illustrated in Fig. 28 by an UrQMD [101] calculation of 
Bass and Dumitru [102] for central Au-I-Au collisions at top RHIC energy. 
We select here the results concerning the survival of the hadronic multiplici- 
ties Ni throughout the dynamics of the hadronic expansion phase, which we 
have postulated above, based on the equality of the hadronization temper- 
atures, Th ~ 160 MeV, observed in e+e^ annihilation (Fig. 19), where no 
hadronic expansion phase exists, and in central collisions oi A m 200 nuclei 
(Figs. 26, 27). In fact. Fig. 28 shows that the {Ni} observed at the end of 
the hadronic cascade evolution agree, closely, with the initial {Ni} as derived 
from a Cooper- Frye procedure (equation 110) directly at hadronization. On 
the other hand, pt spectra and radial flow observables change, drastically, 
during the hadronic cascade expansion phase. 
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Fig. 28. Modification of mid-rapidity hadron multiplicities in central Au+Au colli- 
sions at ^/s — 200 GeV after chemical freeze-out at T = Tc. Squares show a hydro- 
dynamic model prediction at T = Tc (without further interaction); circles show the 
result of an attached UrQMD hadronic cascade expansion calculation [102]. 



The hadronic multiplicity distribution {iV^}, arising from the hadronization 
process at high y/s, freezes-out instantaneously also in A-l-A collisions, and 
is thus preserved throughout the (isentropic) hadronic expansion phase. It is 
thus directly measurable and, moreover, its hadrochcmical equilibrium features 
lend themselves to an analysis within the framework of Hagedorn-type sta- 
tistical, grand canonical models. As we shall show below, the outcome of this 
analysis is contained in a [Tp, /is] parameter pair that reflects the conditions 
of QCD matter prevailing at hadronization, at each considered y^. In fact, 
the [T, /i] points resulting from the SHM analysis exhibited in Figs. 26, 27 
(at ^/s = 17.3 and 200 GeV, respectively) have been shown in the QCD mat- 
ter phase diagram of Fig. 1 to approach, closely, the parton-hadron phase 
coexistence line predicted by lattice QCD. Thus, Th ~ Tc at high y/s: hadro- 
chemical freeze-out occurs in the immediate vicinity of QCD hadronization, 
thus providing for a location of the QCD phase boundary. 



4.2 Grand canonical strangeness enhancement 

The statistical model analysis [21, 95, 96] of the hadronization species distri- 
bution Ni in A-l-A collisions is based on the grand canonical partition function 
for species i, 

iz=^r ^ nn^ 

67T^T Jo E,{k)exp{iE,ik)~,,,)/T}±l ^ ^ 

where Ef — k'^ + mf, and Hi = jiBBi -t- UsSi + njl^ is the total chemical 
potential for baryon number B, strangeness S and isospin 3-component I3. 
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Its role in equation 113 is to enforce, on average over the entire hadron source 
volume, the conservation of these quantum numbers. In fact, making use of 
overall strangeness neutrality {J2i ^iSi = 0) as well as of conserved baryon 
number (participant Z+N) and isospin (participant (N-Z)/Z) one can reduce 
Hi to a single effective potential jib- Hadronic freeze-out is thus captured in 
three parameters, T, V and /Xfe. The density of hadron/resonance species i then 
results as 

T S 

ni = ——lnZi (114) 



V6^ 



which gives 



We see that the common freeze-out volume parameter is canceled if one 

considers hadron multiplicity ratios, Ni/Nj, as was done in Fig. 27. Integration 
over momentum yields the one-particle function 

VTgi 2 (±l)"+i /nmA /n/zA 
= ^ E ^2 (— j exp (— j (116) 

n— 1 

where K2 is the modified Bessel hmction. At high T the; c;ffccts of Bose or 
Fermi statistics (represented by the ±1 term in the denominators of equa- 
tion 113 and equation 115) may be ignored , finally leading to the Boltzmann 
approximation 

VTgi 2 



Ni 



i^.(f)e.p(f) (117) 



2^2 

which is the first term of equation 116. This approximation is employed 
throughout the SHM analysis. It describes the primary yield of hadron species 
i, directly at hadronization. The abundance of hadronic resonance states is ob- 
tained convoluting equation 115 with a relativistic Brcit-Wigner distribution 
[21]. Finally, the overall multiplicity, to be compared to the data, is determined 
as the sum of the primary multiplicity (cf. equation 117) and the contributions 
arising from the unresolved decay of heavier hadrons and resonances: 

^observed ^ ^primary ^ ^ q^^- ^ ^^^g^ 

3 

After having exposed the formal gear of grand canonical ensemble analysis 
we note that equation 117 permits a simple, first orientation concerning the 
relation of T to /i^ in A-l-A collisions by considering, e.g., the antiproton to 
proton production ratio. From equation 117 we infer the simple expression 

N{p)/N{p) = exp{-2fiB/T). (119) 



Taking the mid-rapidity value 0.8 for p/p (from Fig. 27) at top RHIC energy, 
and assuming that hadronization occurs directly at the QCD phase boundary, 
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and hence T ~ 165 MeV, we get fig ~ 18 MeV from equation 119, in 
close agreement with the result, hb = 20 MeV, obtained [96] from the full 
SHM analysis. Equation 119 illustrates the role played by hb in the grand 
canonical ensemble. It logarithmically depends on the ratio of newly created 
quark-antiquark pairs (the latter represented by the p yield), to the total 
number of quarks including the net baryon number-carrying valence quarks 
(represented by the p yield). 




Fig. 29. The ratio of total K'^ to total tt^ multiplicity as a function of -^s, in 
central Au+Au and Pb+Pb collisions and in p+p minimum bias collisions [88]. 



The most outstanding property of the hadronic multiplicities observed in cen- 
tral A-f A collisions is the enhancement of all strange hadron species, by fac- 
tors ranging from about 2 to 20, as compared to the corresponding production 
rates in elementary hadron-hadron (and e+e~ annihilation) reactions at the 
same y/s. I.e. the nuclear collision modifies the relative strangeness output by 
a "nuclear modification factor", i?^^ = Npart ■ Npp, which depends 

on -y/fi and Npart and features a hierarchy with regard to the strangeness 
number s — 1,2,3 of the considered species, Rf^i < R^^z < ^^=3- These 
properties are illustrated in Figs. 29 and 30. The former shows the ratio of 
total to positive pion multiplicities in central Au-|-Au/Pb-|-Pb collisions, 
from lower AGS to top RHIC energies, in comparison to corresponding ratios 
from minimum bias p-|-p collisions [88]. We have chosen this ratio, instead 
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of (-ftr+) /Npart, because it reflects, rather directly, the "Wroblewski ratio" 
of produced strange to non-strange quarks [95], contained in the produced 
hadrons, 

{u) + (d) + {u) + (d) [0.45 in AA. ^ ' 

The low value of As in pp (and all other elementary) collisions reflects a quark 
population far away from m, d, s flavor equilibrium, indicating strangeness sup- 
pression [97]. 
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Fig. 30. The nuclear modification factors -Rs=i,2,3 for hyperon and anti-hyperon 
production in nucleus-nucleus collisions at y^s — 17.3 GeV, relative to the p-l-p 
reference at the same energy scaled by Nw{= Npart) - The NA49 data refer to total 
4n yields [103]. Closed lines represent the inverse strangeness suppression factors 
from ref. [104], at this energy. 



The so-called strangeness enhancement property of A-l-A collisions (obvious 
from Figs. 29 and 30) is, thus, seen as the removal of strangeness suppression; 
it is also referred to as a strangeness saturation., in SHM analysis [95, 96], 
for the reason that « 0.45 corresponds to the grand canonical limit of 
strangeness production, implicit in the analysis illustrated in Figs. 26 and 27. 
The average Rf^i at y/s > 10 GeV thus is about 2.2, both in the data of 
Fig. 29 and in the statistical model. It increases (Fig. 30) toward about 10 in 
s = 3 production of f2 hyperons. 
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In order to provide for a first guidance concerning the above facts and termi- 
nology regarding strangcnc!ss production wc propose an extremely naive ar- 
gument, based on the empirical fact of a universal hadronization temperature 
(Figs. 19, 26, 27) at high ^/s. Noting that (s) = (s) and (u) « (u) « (d) « (d) 
in a QGP system at /i5 near zero, and T = 165 MeV, just prior to hadroniza- 
tion, As « (s) / (u) ~ exp{{mu — ms)/T} = 0.45 at ^ if we take current 
quark masses, rus — m„ « 130 MeV. I.e. the value of As in A+A collisions 
at high ^/s resembles that of a grand canonical QGP at /ifc ^ 0, as was in- 
deed shown in a 3 flavor lattice QCD calculation [105] at T « T^. On the 
contrary, a p+p collision features no QGP but a small fireball volume, at 
T Tc, within which local strangeness neutrality, (s) = (s) has to be strictly 
enforced, implying a canonical treatment [97]. In our naive model the expo- 
nential penalty factor thus contains twice the strangeness quark mass in the 
exponent. As in pp collisions w exp{2{mu — ms)/T} w 0.2, in agreement with 
the observations concerning strangeness suppression, which are thus referred 
to as canonical suppression. In a further extension of our toy model, now ig- 
noring the u, d masses in comparison to TOs w 135 MeV , we can estimate the 
hierarchy of hyperon enhancement in A-l-A collisions, 

R^^ (X N^'^/N^P-Q.b Npart ~ exp{{-sms + 2sm,)/T} = 2.2, 5.1, 11.6 (121) 

for s = 1, 2, 3, respectively. Fig. 30 shows that these estimates correspond well 
with the data [103] for R^^ derived in An acceptance for A, E and /? as well 
as for their antiparticles, from central Pb+Pb collisions at -y/s = 17.3 GeV. 
The p-hp reference data, and C-l-C, Si-I-Si central collisions (obtained by frag- 
mentation of the SPS Pb beam) refer to separate NA49 measurements at the 
same energy. 

The above, qualitative considerations suggest that the relative strangeness 
yields reflect a transition concerning the fireball volume (that is formed in 
the course of a preceding dynamical evolution) once it enters hadronization. 
Within the small volumes, featured by elementary collisions (see section 4.3), 
phase space is severely reduced by the requirement of local quantum number 
conservation [97, 106] including, in particular, local strangeness neutrality. 
These constraints are seen to be removed in A-|-A collisions, in which ex- 
tended volumes of high primordial energy density are formed. Entering the 
hadronization stage, after an evolution of expansive cooling, these extended 
volumes will decay to hadrons under conditions of global quantum mechanical 
coherence, resulting in quantum number conservation occuring, non-locally, 
and on average over the entire decaying volume. This large coherent volume 
decay mode removes the restrictions, implied by local quantum number bal- 
ancing. In the above, naive model we have thus assumed that the hadroniza- 
tion of an Omega hyperon in A+A collisions faces the phase space penalty 
factor of only three .s quarks to be gathered, the corresponding three s quarks 
being taken care of elsewhere in the extended volume by global strangeness 
conservation. In the framework of the SHM this situation is represented by 
the grand canonical ensemble (equations 115, 117); the global chemical po- 
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tential expresses quantum number conservation on average. Strict, local 
conservation is represented by the canonical ensemble. 

The grand canonical (GC) situation can be shown to be the large collision 
volume limit (with high multiplicities {Ni}) of the canonical (C) formulation 
[107, 104], with a continuous transition concerning the degree of canonical 
strangeness suppression [104]. To see this one starts from a system that is al- 
ready in the GC limit with respect to baryon number and charge conservation 
whereas strangeness is treated canonically. Restricting to s = 1 and — 1 the 
GC strange particle densities can be written (from equation 117) as 

nfiii = Af (122) 

with 

Zs=±i = ^mlK2 (^) exp{iBsfiB + ^Mq)/?^} (123) 

and a "fugacity factor" A^^ = exp {jis/T). The canonical strange particle 
density can be written as [104] 

nf = nf^.(A;) (124) 
with an effective fugacity factor 

where S±i = X]s=±i ■Z^s=±i is the sum over all created hadrons and resonances 
with s = ±1, the In{x) are modified Bessel functions, and x = 2y/ SiS-i is 
proportional to the total fireball volume V. In the limit x w F — > oo the 
suppression factor Ii{x) / Iq{x) 1, and the ratio S±i/ ^ Si S-i corresponds 
exactly to the fugacity Xg in the GC formulation (see equation 122). Thus 
the C and GC formulations are equivalent in this limit, and the canonical 
strangeness suppression effect disappears. Upon generalization to the complete 
strange hadron spectrum, with s = ±1,±2,±3, the strangeness suppression 
factor results [104] as 

rj{s) = h{x)/Io{x). (126) 

In particular for small x (volume), 77(5) —> (x/2)*, and one expects that 
the larger the strangeness content of the particle the smaller the suppression 
factor, and hence the larger the enhancement in going from elementary to cen- 
tral A+A collisions. This explains the hierarchy addressed in equation 121, 
and apparent from the data shown in Fig. 30. In fact, the curves shown in this 
figure represent the results obtained from equation 126, for s = 1, 2, 3 hyperon 
production at ^/s = 17.3 GeV [104]. They are seen to be in qualitative agree- 
ment with the data. However the scarcity of data, existing at top SPS energy 
for total hyperon yields, obtained in 47r acceptance (recall the arguments in 
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section 4.1) both for A+A and p+p collisions does not yet permit to cover 
the SHM strangeness saturation curves in detail, for s > 1. 

This saturation is seen in Fig. 30, to set in already at modest system sizes, 
but sequentially so, for ascending hyperon strangeness. Note that SHM satu- 
ration is sequentially approached, from equation 126, with increasing fireball 
volume V . In order to make contact to the experimental size scaling with cen- 
trality, e.g. Npart, the model of ref. [104], which is illustrated in Fig. 30, has 
converted the genuine volume scale to the Npart scale by assuming a universal 
eigenvolume of 7 fm^ per participant nucleon. I.e. Npart — 10 really means 
a coherent fireball volume of 70 fm^, in Fig. 30. Within this definition, satu- 
ration of s = 1, 2, 3 sets in at fireball volumes at hadronization of about 60, 
240 and 600 fm^, respectively: this is the real message of the SHM curves in 
Fig. 30. 




100 200 300 400 
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Fig. 31. The mid-rapidity if"*" to tt"'" ratio vs. Npart in minimum bias Au-|-Au 
collisions at y/s = 200 GeV, compared to the percolation model [108] (solid line); 
a prediction of which for Cu4-Cu at similar energy is given by the long dashed line 
(see text for detail). 



The above direct translation of coherent fireball volume to participant number 
is problematic [108] as it assumes that all participating nucleons enter into a 
single primordially coherent fireball. This is, however, not the case [108] par- 
ticularly in the relative small scattering systems that cover the initial, steep 
increase of 7y(s), where several local high density clusters are formed, each 
containing a fraction of Npart- This is revealed by a percolation model [108] of 
cluster overlapp attached to a Glauber calculation of the collision/energy den- 
sity. At each Npart an average cluster volume distribution results which can 
be transformed by equation 126 to an average {ri{s, V)} distribution whose 
weighted mean is the appropriate effective canonical suppression factor corre- 
sponding to Npart- On the latter scale, the SHM suppression curve thus shifts 
to higher Npart, as is shown in Fig. 31 for the K'^/n^ ratio vs. Npart, measured 
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at mid-rapidity by PHENIX in Au-I-Au collisions at i/i = 200 GeV, which is 

reproduced by the percolation model [108]. Also included is a prediction for 
Cu+Cu at this energy which rises more steeply on the common Npart scale 
because the collision and energy density reached in central Cu-I-Cu collisions, 
at Npart ~ 100, exceeds that in peripheral Au+Au collisions (at the same 
Npart) which share a more prominent contribution from the dilute surface re- 
gions of the nuclear density profile. We note, finally, that this deviation from 
universal Npart scaling does not contradict the observations of a perfect such 
scaling as far as overall charged particle multiplicity densities are concerned 
(recall Fig. 14) which are dominated by pions, not subject to size dependent 
canonical suppression. 

4.3 Origin of hadro- chemical equilibrium 

The statistical hadronization model (SHM) is not a model of the QCD con- 
finement process leading to hadrons, which occurs once the dynamical cooling 
evolution of the system arrives at Tc- At this stage the partonic reaction 
volume, small in elementary collisions but extended in A-|-A collisions, will 
decay (by whatever elementary QCD process) to on-shell hadrons and res- 
onances. This coherent quantum mechanical decay results in a de-colierent 
quasi- classical, primordial on-shell hadron-resonance population which, at the 
instant of its formation, lends itself to a quasi-classical Gibbs ensemble de- 
scription. Its detailed modalities (canonical for small decaying systems, grand 
canonical for extended fireballs in A-|-A collisions), and its derived parameters 
[T, /xb] merely recast the conditions, prevailing at hadronization. The success 
of SHM analysis thus implies that the QCD hadronization process ends in 
statistical equilibrium concerning the hadron-resonance species population. 

In order to identify mechanisms in QCD hadronization that introduce the 
hadro-chemical equilibrium we refer to jet hadronization in e^e~ annihilation 
reactions, which we showed in Fig. 19 to be well described by the canonical 
SHM. In di-jet formation at LEP energy, = 92 GeV, we find a charged 
particle multiplicity of about 10 per jet, and we estimate that, likewise, about 
10 primordial partons participate on either side of the back-to-back di-jet [74]. 
There is thus no chance for either a partonic or hadronic, extensive rescat- 
tering toward chemical equilibrium. However, in the jet hadronization models 
developed by Amati and Veneziano [72], Webber [109] and Ellis and Geiger 
[74] the period of QCD DGLAP parton shower evolution (and of perturbative 
QCD, in general) ends with local color neutralization, by formation of spa- 
tial partonic singlet clusters. This QCD "color pre- confinement" [72] process 
reminds of a coalescence mechanism, in which the momenta and the initial 
virtual masses of the individual clustering partons get converted to inter- 
nal, invariant virtual mass of color neutral, spatially extended objects. Their 
mass spectrum [109] extends from about 0.5 to 10 GeV. This cluster mass 
distribution, shown in Fig. 32, represents the first stochastic element in this 
hadronization model. 
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Fig. 32. Invariant mass spectrum of color neutralization clusters in the Veneziano- 
Webber hadronization model [72, 109]. 



The clusters are then re-interpreted within non-perturbative QCD: their 
internal, initially perturbative QCD vacuum energy gets replaced by non- 
perturbative quark and gluon condensates, making the clusters appear like 
hadronic resonances. Their subsequent quantum mechanical decay to on-shell 
hadrons is governed by the phase space weights given by the hadron and 
resonance spectrum [74, 109]. I.e. the clusters decay under "phase space dom- 
inance" [74] , the outcome being a micro-canonical or a canonical hadron and 
resonance ensemble [73, 95]. The apparent hadro-chemical equilibrium thus is 
the consequence of QCD color neutralization to clusters, and their quantum 
mechanical decay under local quantum number conservation and phase space 
weights. We note that the alternative description of hadronization, by string 
decay [110], contains a quantum mechanical tunneling mechanism, leading to 
a similar phase space dominance [111]. 

Hadronization in e~^e~ annihilation thus occurs from local clusters (or 
strings), isolated in vacuum, of different mass but similar energy density cor- 
responding to QCD confinement. These clusters are boosted with respect to 
each other but it was shown [112] that for a Lorentz invariant scalar, such 
as multiplicity, the contributions of each cluster (at similar T) can be repre- 
sented by a single canonical system with volume equal to the sum of clusters. 
In the fit of Fig. 19 this volume sum amounts to about 25 fm^ [73]; the in- 
dividual cluster volumes are thus quite small, of magnitude a few fm^ [74]. 
This implies maximum canonical strangeness s\.ippression but may, in fact, 
require a micro-canonical treatment of strangeness [97], implying a further 
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suppression. These MC effects are oftentimes included [113] in the canonical 
partition functions by an extra strangeness fugacity parameter 7^ < 1 which 
suppresses s = 1,2, 3 in a hierarchical manner, (A^j(s)) « (75)***. The fit of 
Fig. 19 requires 7s = 0.66, a value typical of canonical multiplicity analysis in 
p + p, p + p and e+e" annihilation collisions [97] at ^/s > 30 GeV. 

The above picture, of hadrochemical equilibrium resulting from the com- 
bined stochastic features of QCD color neutralization by cluster formation, 
and subsequent quantum mechanical decay to the on-shcU hadron and reso- 
nance spectrum (under phase space governance) lends itself to a straight for- 
ward extension to A-|-A collisions. The essential new features, of grand canon- 
ical hadronization including strangeness enhancement, should result from the 
fact that extended space-time volumes of e > Ccrit are formed in the course of 
primordial partonic shower evolution, an overlap effect increasing both with 
^/s and with the size of the primordial interaction volume. As the volume of 
the elementary hadronization clusters amounts to several fm^ it is inevitable 
that the clusters coalesce, to form extended "super-cluster" volumes prior to 
hadronization [108]. As these super-clusters develop toward hadronization via 
non perturbative QCD dynamics, it is plausible to assume an overall quantum 
mechanical coherence to arise over the entire extended volume, which will 
thus decay to hadrons under global quantum number conservation, the decay 
products thus modeled by the GC ensemble. 

Our expectation that space-time coalescence of individual hadronization 
clusters will lead to a global, quantum mechanically coherent extended super- 
cluster volume, that decays under phase space dominance, appears as an anal- 
ogy to the dynamics and quantum mechanics governing low energy nuclear 
fission from a preceding "compound nucleus" [114]. Note that the observation 
of a smooth transition from canonical strangeness suppression to grand canon- 
ical saturation (Figs. 30, 31) lends further support to the above picture of a 
percolative growth [108] of the volume that is about to undergo hadronization. 

An extended, coherent quark gluon plasma state would, of course, repre- 
sent an ideal example of such a volume [115] and, in fact, we could imagine 
that the spatial extension of the plasma state results from a percolative over- 
lap of primordial zones of high energy density, which becomes more promi- 
nent with increasing ^/s and Npart- A QGP state preceding hadronization will 
thus lead to all the observed features. However, to be precise: the hadroniz- 
ing QCD system of extended matter decaying quantum coherently, could still 
be a non-equilibrium precursor of the ideal equilibrium QGP, because we 
have seen above that hadrochemical equilibrium also occurs in e~^e~ annihi- 
lation, where no partonic equilibrium exists. It gets established in the course 
of hadronization, irrespective of the degree of equilibrium prevailing in the 
preceding partonic phase. 
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4.4 Hadronization vs. rapidity and y/s 

We have argued in section 4.1 that, at relatively low ^/s, the total rapidity 
gap Ay does not significantly exceed the natural thermal rapidity spreading 
width Fi ss 2.35 {r/m.i)'^/'^ of a single, isotropically decaying fireball, cen- 
tered at mid-rapidity and emitting hadrons of mass mi [98]. However, this 
procedure involves an idealization because in the real Pb+Pb collision the 
intersecting dilute surface sections of the nuclear density profiles will lead 
to a significant contribution of single-scattering NN collisions, outside the 
central high density fireball. The leading hadron properties of such "corona 
collisions" result in wider proper rapidity distributions, quite different from 
those of the central fireball decay hadrons. Their contribution will thus be 
prominent near target/projectile rapidity, and will feature a canonically sup- 
pressed strangeness. The one-fireball assumption, although inevitable at small 
Ay, does not quite resemble the physical reality. This may explain the need 
for an additional strangeness suppression factor in the GC one-particle par- 
tition function (equation 113) that has, unfortunately, also been labeled 7., 
but expresses physics reasons quite different from the extra suppression factor 
that reflects micro-canonical phase space constraints in elementary collisions. 
It turns out that all GC analysis of central A+A collisions at low y^, and 
addressed to total Att multiplicities, requires a 7^ of 0.7 - 0.85 [21]; in the fit 
of Fig. 26 7« = 0.84. 

At RHIC, Ay « 11 ^ Fi, and such difficulties disappear: 73 « 1 at midra- 
pidity and, moreover, the wide gap permits a SHM analysis which is differen- 
tial in y. Fig. 33 shows the y-dependence of the ratios 7r~/7r+, K~ /K^ andp/p 
as obtained by BRAHMS [116] in central Au+Au coUisions at = 200 GeV. 
The figure shows a dramatic dependence of the p/p ratio, which reflects the 
local baryochemical potential according to equation 119. At ycM > 1 the p/p 
ratio drops down steeply, to about 0.2 at y « 3.5, thus making close con- 
tact to the top SPS energy value obtained by NA49 [117]. The K~/K+ ratio 
follows a similar but weaker drop-off pattern, to about 0.65 again matching 
with the top SPS energy value of about 0.6 [118]. The deviation from unity of 
these ratios reflects the rapidity densities of initial valence u,d quarks, rela- 
tive to the densities of newly created light and strange quark-antiquark pairs, 
i.e. the y distribution of the net baryon number density, and of the related 
baryo-chemical potential of the GC ensemble. Thus, in analyzing successive 
bins of the rapidity distributions in Fig. 33, the major variation in the GC 
fit concerns the baryo-chemical potential /i_B(y) which increases from about 
20 MeV (Fig. 27) at mid-rapidity, to about 150 MeV at y > 3 while the 
hadronization temperature stays constant, at T = 160 MeV. This interplay 
between K^/K'^, p/p and is illustrated [116] in the right hand panel of 
Fig. 33, and shown to be well accounted for by the GC statistical model [119]. 

These considerations imply that hadronization at RHIC (and LHC) energy 
occurs local in j/-space and late in time. The density distribution of net baryon 
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Fig. 33. (left) Anti-hadron to hadron ratios as a function of rapidity in central 
Au+Au collisions at -^s = 200 GeV. The right panel shows an interpretation of 
the correlation between p/p and K~ /K^ in terms of baryo-chemical potential fj,B 
variation in the grand canonical statistical model. From [116]. 



number results from the primordial pQCD shower evolution (c.f. section 3.4), 
and is thus fixed at formation time, to < 0.5 fm/c at RHIC. Hadronization 
of the bulk partonic matter occurs later, at t > 3 fm/c [75, 84], and trans- 
mits the local conditions in rapidity space by preserving the local net baryon 
quantum number density. Most importantly we conclude that hadronization 
occurs, not from a single longitudinally boosted fireball but from a succession 
of " super-clusters" , of different partonic composition depending on y, and de- 
caying at different time due to the Lorentz-boost that increases with y, in an 
"inside-outside" pattern (c.f. Fig. 2). We are thus whitnessing at hadroniza- 
tion a Hubble expanding system of local fireballs. The detailed implications 
of this picture have not been analyzed yet. Note that a central RHIC collision 
thus does not correspond to a single hadronization " point" in the [T, /i] plane 
of Fig. 1 but samples {T, /x} along the QCD parton-hadron coexistence line 
[120]. 

Throughout this section we have discussed hadronic freeze-out at high 
a/s only (top SPS to RHIC energy), because of the proximity of the chem- 
ical freeze-out parameters [T, iih] to the QCD phase boundary from lattice 
QCD, which suggests an overall picture of hadronization, to occur directly 
from a partonic cluster or super-cluster. Our discussion of the GC statistical 
hadronization model has been explicitly or implicitly based on the assump- 
tion that hadronic freeze-out coincides with hadronization. However, the GC 
model has also been applied successfully to hadro-chemical freeze-out at yfs 
down to a few GeV [21, 95, 96] where it is not expected that the dynamical 
evolution traverses the phase boundary at all, but grand canonical multiplic- 
ity distributions, and their characteristic strangeness enhancement pattern, 
are observed throughout. Toward lower ^/s, T decreases while increases, 
as is shown in Fig. 34 which presents a compilation of all reported freeze-out 
parameters [96]. 
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Fig. 34. Energy dependence of the hadro-chemical freeze-out points obtained by 
grand canonical statistical model analysis in the plane [T, ^_g], with interpolating 
curve at fixed energy per particle of about 1 GeV [95, 121]. 



These points have also been included in the phase diagram of Fig. 1 which 
shows that they are gradually branching away from the phase separation 
boundary Hne that could recently be predicted by lattice QCD once new 
methods had been developed to extend the theory to finite jiB [H, 12]. At 
y/s > 20 GeV we see that 

e,{QCD) ^CH- ecc (127) 

where ecc is the freeze-out density inferred from GC analysis [21, 95, 96]. 

In turn, the GC hadronic freeze-out points drop below the lattice QCD 
coexistence line at lower ^/s^ implying that chemical freeze-out now occurs 
within the hadronic expansion phase. This requires a model of freeze-out, now 
governed by the properties of a high density hadronic medium, upon expansive 
cooling and dilution. Holding on to the model of a quantum mechanical de- 
coherence decay to on-shell hadrons that we discussed in section 4.3, we argue 
that an initial, extended high density hadronic fireball, given sufficient life- 
time at T smaller, but not far below Tcj could also be seen as a quantum 
mechanically coherent super-cluster, as governed by effective mean fields [122]. 
In such a medium hadrons, at T near X'^, acquire effective masses and/or 
decay widths far off their corresponding properties in vacuum: they are off- 
shell, approaching conditions of QCD chiral symmetry restoration as T —> Tc 
[123]. This symmetry is inherent in the elementary QCD Langrangian, and 



Global Properties of Nucleus-Nucleus Collisions 79 

"softly" broken within the light quark sector by the small non-zero current 

quark masses, but severely broken at T ^ by the high effective constituent 
quark masses that get dressed by non perturbative QCD vacuum condensates. 
Pictorially speaking, hadrons gradually loose this dressing as T —> Tc [124], 
introducing a change, away from in vacuum properties, in the hadronic mass 
and width spectrum. Such in-medium chiral restoration effects have, in fact, 
been observed in relativistic A-l-A collisions, by means of reconstructing the 
in-medium decay of the p vector meson to an observed e^e'^ pair [125]. 

A dense, high T hadronic system, with mean-field induced off-shell con- 
stituents is also, clearly, quantum mechanically coherent. At a certain charac- 
teristic density, e < ec, and temperature T < T,., as reached in the course of 
overall hadronic expansion, this extended medium will undergo a decoherence 
transition to classical on-shell hadrons. Its frozen-out hadronic multiplicity 
distribution should be, again, characterized by the phase space weights of 
a grand canonical ensemble at T < Tg. Theoretical studies of such a mean 
field hadronic expansion mode [126] have also shown that such mechanisms 
play essentially no role at ^/s > 20 GeV because the expanding system is 
already in rapid flow once it traverses the phase boundary, with an expansion 
time scale shorter than the formation time scale of mean field phenomena. At 
lower energies, on the other hand, the system might not even dive into the 
deconfined phase but spend a comparatively long time in its direct vicinity, 
at the turning point between compression and re-expansion where all dynam- 
ical time constants arc large, and the hadron density is high, such that the 
inelastic hadronic transmutation rate becomes high (particularly in collisions 
of more than two hadronic reactants, with reaction rates [127] proportional 
to e"), and sufficiently so for maintaining hadronic chemical equilibrium after 
it is first established at maximum hadron density, in low y/s systems that do 
not cross the phase boundary at all. 

The GC freezc-o\it parameters [T, /i] at various in Fig. 34 permit a 
smooth interpolation in the T,jj plane [121], which, in turn, allows for GC 
model predictions which are continuous in ^/s. Such a curve is shown in Fig. 29 
compared to the An data points for the /n^ multiplicity ratio in central 
collisions Au+Au/Pb+Pb, at all ^/s investigated thus far. It exhibits a smooth 
maximum, due to the interplay of T saturation and /is fall-off to zero, but 
does not account for the sharp peak structure seen in the data at ^/s w 7 GeV 
and jiB ^ 480 MeF. This behavior is not a pecularity of the channel only; 
it also is reflected in an unusually high Wroblewski ratio (see equation 120) 
obtained at = 1 .QGeV , of = 0.60 [21]. This sharp strangeness maximum 
is unexplained as of yet. It implies that hadron formation at this ^ reflects 
influences that are less prominent above and below, and most attempts to 
understand the effect [128, 129, 130] are centered at the assumption that at 
this particular ^ the overall bulk dynamics will settle directly at the phase 
boundary where, moreover, finite jiB lattice theory also expects a QCD critical 
point [11, 12, 13]. This would cause a softest point to occur in the equation of 
state, i.e. a minimum in the relation of expansion pressure vs. energy density, 
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slowing down the dynamical evolution [131, 132], and thus increasing the 

sensitivity to expansion modes characteristic of a first order phase transition 
[130], which occurs at > 1^1"'^ ■ Such conditions may modify the K/tt ratio 
(Fig. 29) [130]. 

It thus appears that the interval from top AGS to lower SPS energy, 
5 < -^i < 10 GeV, promises highly interesting information regarding the 
QCD phase diagram (Fig. 1) in the direct vicinity of the parton-hadron co- 
existence line. In particular, the physics of a critical point of QCD matter 
deserves further study. Observable consequences also comprise so-called "crit- 
ical fluctuations" [133, 134] of multiplicity density, mean transverse momen- 
tum and hadron-chemical composition [135], the latter in fact being observed 
near ^/s = 7 GeV in an event by event study of the Kj-n ratio in central 
Pb-hPb collisions [136]. 
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